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ABSTRACT 


PRC  Information  Sciences  Company  Report  R-1939 
Program  Maintenance  for  the  Reference  Scene  Software 
Carolyn  B.  Shelton,  October  1976  (Unclassified) 


The  U.S.  Army  Engineer  Topographic  Laboratories  (ETL)  at  Fort  Belvoir  is 
presently  engaged  in  a concentrated  effort  aimed  at  developing  a methodology 
for  generating  radar  reference  scenes  from  raw  cartographic  data.  A central 
part  of  this  effort  is  the  identification  of  the  minimum  set  of  radar 
attributes  required  in  such  reference  scenes  in  order  that  they  provide 
sufficient  correlation  when  compared  with  actual  PPI  radar  images.  The 
objective  of  the  effort  described  herein  is  the  development  of  Reference 
Scene  Generation  Software  (RSS)  to  be  used  by  ETL  as  a research  tool  in 

the  development  of  the  final  reference  scene  generation  criteria. 

The  RSS  is  based  on  Digital  Radar  Landmass  Simulation  Software  (DRI/1S) 
provided  by  the  Naval  Training  Equipment  Center  at  Orlando,  Florida.  The 
- first  step  in  the  development  of  RSS  was  to  convert  these  programs  to  run 
in-house  on  ETL' s CDC  6400.  New  input  routines  were  written  to  permit  the 
use  of  in-house  data  bases,  and  a new  output  routine  was  written  to  permit 
the  radar  scene  output  to  be  displayed  on  ETL’s  DICOMED  plotter. 


The  second  phase  of  the  program  involved  further  modifications  of  the  pro- 
grams to  make  them  more  suitable  for  correlation  work.  Among  the  improve- 
ments added  were  the  capabilities  to  vary  image  resolution,  size  and 
coloring.  The  software  was  also  analysed  and  corrected  to  improve  its 
geometric  accuracy. 


Finally,  a routine  was  added  to  permit  the  incorporation  of  the  altitude 
layover  effect  into  the  reference  scenes.  This  effect  produces  a non- 
uniform  radial  translation  of  the  points  on  the  radar  image  and  its 
inclusion  is  expected  to  improve  the  correlation  obtainable  with  the 
reference  scenes.  , ' 
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PREFACE 


This  work  was  principally  performed  by  Ms.  Carolyn  B.  Shelton, 
Planning  Research  Corporation,  under  contract  number  DAAK02-75-C-0098, 
Radar  Programs  Conversion  for  the  U.S.  Army  Engineer  Topographic 
Laboratories,  Fort  Belvoir,  Virginia,  Bruce  B.  Zimmerman,  Contracting 
Officer's  Technical  Representative. 

This  document  is  a Program  Maintenance  Manual  for  the  Reference 
Scene  Software  (RSS)  presently  being  used  by  the  USAETL  in  their 
radar  scene  correlation  studies.  As  such,  it  provides  sufficient 
information  to  enable  the  programmer  to  perform  such  software 
modifications  as  may  be  required  during  the  continuing  development 
of  the  system. 

This  document  should  be  used  in  conjunction  with  the  RSS  User's 
Manual  (ETL-066) • j / 

RSS  is  based  on  Digital  Radar  Landmass  Simulation  Software 
supplied  by  the  Naval  Training  Equipment  Center,  Orlando,  Florida. 
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1. 


General  Description 


1.1  System  Application.  ETL  is  presently  engaged  in  a concentrated 
effort  aimed  at  developing  a methodology  for  generating  radar  refer- 
ence scenes  from  raw  cartographic  data.  A central  part  of  this  effort 
is  the  identification  of  the  minimum  set  of  radar  attributes  required 
in  such  reference  scenes  in  order  that  they  provide  sufficient  corre- 
lation when  compared  with  actual  PPI  radar  images.  RSS  is  being 

used  by  ETL  as  a research  tool  in  the  development  of  such  a set  of 
reference  scene  generation  criteria. 

1.2  Equipment  Environment.  The  software  described  herein  is  written 
in  FORTRAN  IV  for  a CDC  6400  large-scale  computer  using  the  SCOPE  3.4 
operating  system.  Because  of  the  inherently  large  storage  require- 
ments of  the  programs,  on-line  disk  storage  is  used  extensively. 

A DICOMLD  Model  36  scanner/plotter  is  presently  being  used  as  the 
output  display  device. 
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2.  Detailed  Software  Description 
2.1.  System  Description 


The  radar  scene  simulation  process  employed  by  this  system  consists 
of  four  subsystems:  (1)  the  planimetry  preparation  subsystem,  (2)  the 

terrain  preparation  subsystem,  (3)  the  scene  generation  subsystem,  and 
(4)  the  image  processing  subsystem. 

Programs  RSS1  through  RSS5  comprise  the  planimetry  preparation 
subsystem. 

The  input  to  RSS1  is  a raw  planimetry  data  file  which  contains  the 
X-Y  pairs  defining  the  boundary  of  each  feature  and  a feature  code  identify- 
ing the  feature  type  (e.g.,  road,  lake,  etc.).  This  information  is  recorded 
by  the  ETL  Bendix  Programmable  Digitizer.  RSS1  unpacks  the  feature  infor- 
mation, throwing  away  the  digitizer  command  codes,  and  outputs  the  unpacked 
feature  data  onto  a disk  file.  This  disk  file  is  the  input  to  program  RSS2. 

In  RSS2  the  feature  data  is  converted  from  X-Y  pairs  defining  a 
feature  perimeter,  to  horizontal  strip  data.  This  procedure  result  is  a 
file  with  enclosed  feature  data.  Each  horizontal  strip  is  defined  by  the 
X,  Y coordinates  of  its  left  end,  and  by  a delta-X  value  which  is  the 
strip  length. 

Because  of  limited  core  memory  available,  it  is  necessary  to  divide 
the  map  into  smaller  regions  for  processing.  Program  RSS3  does  this  by 
dividing  the  planimetry  strips  generated  by  RSS2  at  fixed  region  boundaries. 

Program  RSS4  is  a sorting  routine  which  reorders  the  planimetry  data 
generated  by  RSS3.  The  output  of  RSS4  is  an  equivalent  data  base  in 
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which  the  data  is  ordered  by  the  following  priorities  respectively: 

1 - region  number 

2 - merge  priority 

3 - Y coordinate 

4 - X coordinate 

Program  RSS5  merges  the  planimetry  data,  output  from  RSS4,  in 
accordance  with  the  above  priority  list.  The  output  of  RSS5  is  a file 
in  which  the  contents  of  each  region  are  described  by  non-overlapping 
horizontal  strips.  In  order  to  accomplish  this,  strips  from  low 
merge  priorities  (such  as  cities)  are  written  first  and  are  written  over 
where  necessary  by  strips  from  features  with  higher  merge  priorities  (such 
as  lakes) . The  completion  of  this  program  marks  the  end  of  the  planimetry 
preparation  subsystem. 

Program  RSS6  is  the  terrain  preparation  program.  The  terrain  data 
is  divided  into  regions  with  the  same  boundaries  as  that  generated  by 
RSS3. 

Programs  RSS7  through  RSS9  comprise  the  scene  generation  subsystem. 

Program  RSS7  merges  the  planimetry  and  terrain  files  for  those 
regions  lying  within  the  radar  ground  range  of  the  target.  The  output 
from  RSS7  consists  of  one  record  for  each  region  within  the  radar  range 
of  the  target.  Each  record  consists  of  an  appropriate  header,  followed 
by  the  planimetry/terrain  content  of  each  point  in  the  region. 

Program  RSS8  uses  the  cartesian  data  base  output  from  RSS7  to 
construct  a series  of  radial  scan  lines  required  to  simulate  a Plan- 
Position-Indicator  (PPI)  radar  scene.  This  operation  consists  of  cartesian- 
to-polar  coordinate  transformation. 
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The  output  from  RSS8  consists  of  N scan  lines  where  360  /N  is  the  angular 
increment  between  scan  lines.  The  single  record  for  each  scan  line  con- 
tains the  planimetry/terrain  data  for  all  points  lying  along  the  scan  line, 
beginning  at  the  target  location  and  ending  20  nautical  miles  out. 

Program  RSS9  applies  the  radar  effects  to  the  radial  scan  line  data 
from  RSS8.  The  effects  considered  are  discussed  later  under  the  program 
description  for  RSS9.  This  program  also  calls  subroutines  which  scale 
the  image  size  and  convert  the  data  back  to  cartesian  coordinates.  Program 
RSS9  outputs  a data  record  for  each  cartesian  coordinate  to  a disk  file. 
RSS9  marks  the  end  of  the  scene  generation  subsystem. 

At  this  point  the  radar  scene  is  complete.  The  image  processing 
subsystem  begins  with  the  SORT  program  which  orders  the  data  records 
generated  by  RSS9  by  Y coordinate  and  X coordinate. 

Program  RSS10  formats  the  final  image  data  onto  an  output  tape  using 
the  cartesian  coordinate  data  from  program  SORT.  This  output  tape  is  then 
read  into  a DICOMED  plotter  which  outputs  the  radar  image. 

Operationally  this  simulation  system  generates  four  scenes  for  each 
target.  Each  scene  depicts  the  radar  image  of  a target  as  viewed  from  a 
different  altitude.  Programs  RSS1  through  RSS8  need  only  be  run  once 
for  any  given  target.  Programs  RSS9  through  RSS10,  however,  need  to  be 
run  once  for  every  scene  generated. 

The  flow  diagram  for  this  system  is  illustrated  in  Figure  II. 1. 

(1.1  in  User's  Manual) 

2.2  Program  Description 

This  section  gives  a detailed  description  of  each  routine  and  sub- 
routine in  this  system.  Attachments  for  each  program  described  include 
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a program  flowchart,  a program  listing,  and  printouts  by  the  program, 
if  any. 

2.2.1 

2. 2. 1.1  Program  Name.  RSS1 

2. 2. 1.2  Storage.  Approximately  43K  words  of  core. 

2.2.1. 3 Run  Time.  Approximately  92  cp  seconds  and  4 minutes  wall  clock 
time. 

2. 2. 1.4  Cost.  Approximately  24  dollars. 

2. 2. 1.5  Program  Function.  RSS1  unpacks  the  ETL  vector  digitized 
feature  information  and  throws  away  the  digitizer  command  codes.  It  also 
performs  certain  preprocessing  functions  to  be  described  later. 

2. 2. 1.6  Program  Description.  For  each  feature,  RSS1  begins  by  pulling 
off  and  reconstructing  the  feature  code  and  (X,Y)  data  pairs.  All  other 
information  on  the  data  tape  is  discarded.  The  X and  Y offsets  are  sub- 
tracted from  the  (X,Y)  pairs,  and  any  point  having  a negative  coordinate 
is  discarded  since  it  lies  off  of  the  adjusted  map  sheet.  This  procedure 
is  necessary  to  bring  the  coordinate  origin  of  the  planimetry  map  into 
coincidence  with  that  of  the  terrain  rap.  Failure  to  do  this  will  cause  a 
systematic  error  in  the  placement  of  the  planimetry  data. 

The  data  for  areal  features,  such  as  lakes  or  cities,  is  then  checked 
to  see  that  a closed  polygon  is  formed,  and  closure  is  performed  if  required. 
Such  a closure  of  all  polygons  is  required  by  RSS2 . Polygons  lying  par- 
tially off  the  map  sheet  are  closed  along  the  sheet  boundary. 

In  order  to  simplify  later  processing,  the  five  digit  feature  codes 
are  mapped  into  the  set  of  integers  (1,  60).  This  permits  easy  formulation 
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of  translation  tables  which  match  feature  codes  with  feature  attributes 
such  as  merge  priority  and  intrinsic  radar  return  intensity. 

Finally,  RSS1  does  special  processing  on  certain  features.  If  a 
feature  is  a power  line  the  towers  are  separated  from  the  cable  and  each 
piece  is  written  out  as  a separate  feature.  A wide  river  is  represented 
by  two  separate  features,  one  defining  its  right  bank  and  another  defining 
its  left  bank.  In  addition  to  maintaining  the  two  banks,  RSSl  combines 
them  into  a third  feature  which  is  a closed  polygon  representing  the  water 
surface . 

The  following  is  a description  of  the  important  variables  in  RLMS1 . 


Label 

Type 

Description 

Blank 

Integer 

Contains  BCD  representation  for  a 
blank  character. 

Block 

Integer 

Block  counter. 

Buffer 

Integer 

Array 

Digitizer  input  tape  200-word  feature 
record 

Down 

Integer 

Contains  BCD  representation  for  the 
character  "D"  which  represents  the 
down  pen  command  (start  of  a line). 

Etlif 

Integer 

Array 

Array  containing  the  ETL  feature 
codes 

Finish 

Integer 

Contains  BCD  representation  for  the 
character  "B"  wnich  indicates  the 
end  of  clock  . 

Header 

Integer 

Header  information  indicator.  Con- 
tains either  the  BCD  representation 
for  "Yes”  or  "No". 

IByte 

Integer 

Byte  pointer  into  input  buffer. 

I Close 

Integer 

Closed  feature  indicator.  Contains 
either  the  BCD  representation  for 
"Yes"  or  "No". 

IEXP 

Integer 

Exponent  of  the  base  10  required  to 
convert  BCD  character  to  numeric 
machine  code. 

IF 

Integer 

Feature  code 

I MASK 

Integer 

Six  bit  mask 

IN 

Integer 

Word  counter  for  output  buffer 
record 

INDIC 

Integer 

Output  unit  device  check  word 

INDI 

Integer 

Input  unit  device  check  word 

INP1 

Integer 

Value  of  1N+1 

INP2 

Integer 

Value  of  1N+2 

INP3 

Integer 

Value  of  1N+3 

INTAP 

Integer 

Input  unit  device  number.  INTAP=7 

IN  3 

Integer 

Number  of  coordinates  (points  *2)  on 
right  river  bank 

IN3P3 

Integer 

Value  of  IN3+3 

IN4 

Integer 

Number  of  coordinates  (points  *2) 
on  left  river  bank 

INP4P3 

Integer 

Value  of  IN4+3 

ISEQ 

Integer 

Feature  counter 

ISHIFT 

Integer 

Number  of  bits  to  shift  word 

IT 

Integer 

Value  of  IN4+IN3+2  which  equal  the 
total  number  of  coordinates  (points 
in  the  left  and  right  river  banks 

ITEMP 

Integer 

Tinporary  storage  word 

ITYPE 

Integer 

BCD  representation  for  current 
feature  type 

IWD 

Integer 

Number  of  words  in  input  record. 
Initialized  at  180  words 
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KFLAG 

Integer 

Riv'ur  bank  processing  flag. 
1=  one  bank  read 
2=both  banks  read 

KMAX 

Integer 

Total  number  of  words  in  output  record 

KOUNT 

Integer 

Byte  counter  into  next  block 

LOC 

Integer 

Word  pointer  for  input  buffer 

NBYTE 

Integer 

Byte  pointer  for  input  buffer 

NEWBLK 

Integer 

New  block  indicator.  Contains  BCD 
representation  for  either  "Yes"  or  "No" 

NO 

Integer 

Contains  BCD  representation  for  "No" 

NOCLOSE 

Integex 

Array 

Contains  array  of  feature  codes  for 
unclosed  features 

NWD 

Integer 

Word  counter  for  output  buffer 

NWORD 

Integer 

Word  counter  for  input  buffer 

OUTAP 

Ingeger 

Output  unit  device  number.  0UTAP=8 

PEN 

Integer 

Pen  command 

POINT 

Integer 

BCD  representation  for  a decimal  point 

RESTART 

Integer 

Restart  indicator  for  picking  up  the 
remaining  X or  Y coordinate  (BCD  value) 
in  the  next  block.  Contains  BCD  repre- 
sentation for  either  "Yes"  or  "No" 

START 

Integer 

Contains  BCD  representation  for  " ( " 
which  indicates  start  of  new  block 

TEMP 

Real 

Temporary  area  for  extracted  X and  Y 
coordinate  character  values 

TYPE 

Integer 

Array 

Array  containing  BCD  representation 
for  English  description  of  the  feature 
type 

UP 

Integer 

Contains  BCD  representation  for  the 
character  "U"  which  represents  the 
UP  pen  conmand  (end  of  a line) 

rw 
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XMAX 

Real 

Maximum  X coordinate  value  for  this 
feature 

XMIN 

Real 

Minimum  X coordinate  value  for  this 
feature 

X coordinate  offset  (correction  factor) 

XORGIN 

Real 

X2 

Real 

Array 

Feature  output  buffer 

X3 

Real 

Array 

Tenporary  storage  for  right  river 
bank  data 

X4 

Real 

Array 

Temporary  storage  for  left  river  . 

bank  data 

YES 

Integer 

Contains  BCD  representation  for  "No" 

YMAX 

Real 

Maximum  Y coordinate  value  for  this 
feature 

| 

YMIN 

Real 

Minimum  Y coordinate  value  for  this 
feature 

YORGIN 

Real 

Y coordinate  offset  (correction  factor) 

2. 2.1.7  Input . RSS1  requires  one  tape  input  file  (tape  7)  and  one  card 
input  file.  The  tape  input  consists  of  a digitizer  tape  containing  the 
(X,Y)  coordinates  of  the  points  defining  the  location  of  the  planimetry 
features  on  the  map.  The  format  of  this  tape  is  presented  in  Figure  II. 1 
Each  small  box  represents  a 6-bit  BCD  character,  and  each  physical 
record  contains  1800  such  characters  (or  180  CDC  words) . Tbe  beginning 
of  the  information  for  each  feature  is  marked  by  a "(".  Following  this 
are  twelve  characters  which  specify  the  starting  location  for  the  feature 
but  are  not  used.  The  next  five  characters  contain  the  feature  code. 

A tabulation  of  the  possible  feature  codes  is  given  in  Table  II. 1. 
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FEATURE  TYPE  FEATURE  CODE 


RADAR  GREY  SHADES 

20101  THRU  20320 

LARGE  RIVERS  (WATER  ON  LEFT) 

10110 

LARGE  RIVERS  (WATER  ON  RIGHT! 

10123 

OAKS 

10130 

MARSHES  AND  SWAMPS 

10140 

LAKES 

10150 

ISLANDS 

10160 

RIVERS  ANO  STREAMS 

10170 

RAILROAD  YARDS 

10210 

RAILROADS 

10223 

TOWNS  AND  SUBURBS 

10310 

MEDIUM  CITIES  ANO  COMMERCIAL  AREAS 

10320 

BIG  CITIES  ANO  INDUSTRIAL  AREAS 

10330 

LARGE  ISOLATED  BUILDINGS 

10340 

INTERSTATE  HIGHWAYS  AND  TURNPIKES 

10410 

MAJOR  ROADS 

10420 

SECONDARY  ROADS 

10430 

UNPAVED  ROADS  ANO  TRAILS 

10440 

AIRPORT 

10450 

POWER  LINE  TOWERS  (WITH  CABLES! 

10510 

DRIVE-IN  MOVIES 

10520 

FIRE  OR  RADIO  TOWERS 

10530 

CEMETERIES 

10540 

POL  AREA 

10550 

HARDWOOD  FOREST 

10610 

EVERGREEN  FOREST 

10620 

MEADOWS  AND  GRASSY  FIELOS 

10630 

ORY  ROCKY  AREAS 

10640 

SAND  ANO  SAGEBRUSH  AREAS 

10650 

SNOW  COVERED  AREAS 

10660 

ORY  RIVERBEDS*  CANALS,  AND  STORM  ORAINS 

10670 

DRY  LAKE  BEOS  ANO  GULCHES 

10680 

tadle  ii.i  - correspondi:nci:  hktwi:i;n  feature  types 
AND  feature" CODES 

2-11 


Following  the  feature  code  is  a pen-down  command  and  this  signals  the 
beginning  of  the  (X,Y)  data. 

The  (X, Y)  data  for  the  feature  begins  in  the  twentieth  character 
of  the  feature  record.  Each  coordinate  is  in  an  F7.3  format.  This  means 
that  each  digit  (as  well  as  the  decimal  point)  is  represented  as  a 
separate  character,  with  the  maximum  coordinate  value  being  999.999. 

All  coordinates  are  given  in  inches  relative  to  the  table  origin.  The 
last  Y-value  for  the  feature  is  followed  by  a "U"  which  is  the  pen-up 
command  from  the  digitizer. 

A feature  record  may  begin  anywhere  in  a block  of  data,  and  may 
require  more  them  one  block.  Similarly,  a single  block  may  contain  the 
information  for  several  features  of  small  size.  A feature  record  may  be 
broken  at  any  point  except  in  the  middle  of  the  string  of  characters 
representing  the  feature  code  or  a coordinate.  A "B"  is  used  to  signal 
the  end-of-information  for  a particular  block. 

The  card  input  to  RSS1  consists  of  a single  card  containing  the 
X and  Y offsets.  These  numbers  are  subtracted  from  every  (X,Y)  pair 
for  every  feature  and  are  required  to  bring  the  origin  of  the  planimetry 
data  into  coincidence  with  the  origin  used  for  the  terrain  data.  The 
information  is  entered  in  a 2F10.3  format. 

2. 2. 1.8  Output.  RSS1  outputs  the  processed  planimetry  file  to  disk. 
Each  feature  is  represented  by  a record  with  the  following  format: 

Word  1:  Feature  code 

Word  2:  Total  number  of  coordinates  in  the  record  i.e.  twice 

the  number  of  points  N. 
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Word  3-word  2N+2:  X,Y  data 

RLMS1  also  outputs  a printout  identifying  the  type  and  location 
of  all  features  on  the  map  sheet.  A sample  of  this  printout  is  given  in 
Figure  II. 1. 

2. 2. 1.9  Externals.  System  routines  called  by  RSS1  are  LENGTH  and  UNIT 
which  are  mass  storage  I/O  function  routines.  There  are  no  subroutine 
calls  from  this  program. 

2.2.1.10  Error  Conditions.  There  are  two  error  conditions  for  RSS1. 

One  error  condition  occurs  when  there  is  a format  error  found  on  the 

c^9itizer  tape.  The  messages  output  for  this  error  condition  are  as 
follows : 

INCOMPLETE  FEATURE  HEADER  IN  BLOCK  XXXX 
RUNSTREAM  ANALYSIS  TERMINATED 
OR 

ILLEGAL  PEN  COMMAND  IN  BLOCK  XXX 
RUNSTREAM  ANALYSIS  TERMINATED 

The  octal  data  for  the  block  containing  the  error  is  printed  after 

each  of  the  above  two  error  messages. 

The  other  error  condition  occurs  when  an  error  is  found  while 

reading  or  writing  to  or  from  the  input  tape  or  the  output  disk  storage. 

The  messages  generated  for  this  error  condition  are  as  follows : 

I/O  ERROR  DETECTED  AT  BLOCK  XXXX  TOTAL  WORIS 
TRANSFERRED  = XXX  BUFFER  STATUS  = XX 
OR 

OUTPUT  ERROR  DETECTED  AT  BLOCK  XXXX  TOTAL 
WORDS  OUTPUT  = XXX  BUFFER  STATUS  = XX 

2.2.1.11  Program  Flowchart.  The  flow  diagram  for  RSS1  is  illustrated 
in  Figure  II.  2. 

2.2.1.12  Program  Listing.  The  program  listing  for  RSS1  is  attached. 


I 
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FIGURE  II. 2 - PROGRAM  RSS1  FLOWCHART 
(Page  4 of  4) 
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2.2.2 

2. 2. 2.1  Program  Name.  RSS2 

2. 2. 2. 2 Storage.  Approximately  48K  words  of  core. 

2. 2. 2. 3 Run  Time.  Approximately  117  cp  seconds  and  50  minutes  wall 
clock  time. 

2. 2. 2. 4 Cost.  Approximately  350  dollars. 

2. 2. 2. 5 Program  Function.  RSS2  converts  the  X-Y  pairs  describing  each 
feature  boundary  into  horizontal  strip  format.  This  is  done  for  both  "line" 
features  such  as  roads  and  "area"  features  such  as  lakes.  Clearly  a point 
feature  such  as  a building  can  be  displayed  by  a single  strip  of  unit 
length.  For  area  features,  long  horizontal  "fill"  strips  are  generated 

so  that  the  area  enclosed  by  the  feature  boundary  is  now  part  of  the 
feature. 

2. 2. 2.6  Program  Description.  The  process  of  strip  generation  takes  place 

in  three  different  stages.  These  stages  are  point  processing,  line  pro- 
cessing, and  polygon  processing.  Each  feature  is  defined  as  either  a 
point,  line,  or  polygon  feature.  If  the  feature  is  a point  feature,  then 
a six  word  unit  strip  is  formatted  for  each  point  of  the  feature  into  am 
output  buffer  in  the  order  that  they  are  input  by  RSS1.  These  unit  strips 
are  then  ordered  by  Y and  then  X.  After  sorting  is  complete,  a record 
header  is  added  to  the  buffer  (array  IX)  and  the  buffer  is  output  to  a 
disk  file.  For  line  processing,  the  following  procedure  takes  place.  For 
a consecutive  point  pair,  say  and  (X2>Y2),  Ax  and  Ay  is  confuted. 

To  construct  the  data  strips  corresponding  to  these  two  points  begin  by 
calculating  the  inverse  slope. 

s-Ax  - I X.-X/Y  -Y.  I 
£Y 
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This  slope  is  then  scaled  to  a 2048  grid.  S 


2048 

nr' 


If  £Y  is  not  9,  then  S becomes:  S = S +(2*2048) 

Note:  + depends  on  direction  of  slope 

The  number  of  strips  to  be  generated  is  given  by  f Y^-Y^  ' eac^ 

strip  being  of  unit  width. 

For  the  first  strip,  the  values  computed  are  Y=Y , x=X,  and  Ax=X+S . 
For  each  of  the  remaining  strips,  Y=Y^-n,  X=X1+S+n,  and  Ax=S  where  n equals 

the  number  of  strips  processed  for  this  point  pair.  The  Y,  X,  and  Ax 
values  for  each  strip  is  formatted  into  a work  area  (array  JX) . After  the 
strips  are  generated  for  this  point  pair,  the  recently  generated  strips  in 
the  JX  array  are  sorted  into  the  previously  generated  strips  of  the  IX 
array  ordered  by  Y,  X, Ax,  9 respectively.  The  above  procedure  is  repeated 
for  each  point  pair  of  the  feature.  Polygon  processing  follows  line  pro- 
cessing. If  the  line  feature  already  processed  is  determined  to  be  a non- 
polygon or  filled  polygon  with  no  overlaps,  the  feature  record  with  its 
header  added  is  output  to  a disk  file.  Otherwise  a series  of  tests  are 
made  for  each  pair  of  line  segments  (strips)  in  the  feature  to  determine 
where  gaps  or  overlaps  occur.  If  a gap  exists  between  a pair  of  line 
segments,  a filler  strip  is  formatted  and  inserted  in  the  proper  order 
into  the  output  buffer  (IX  array) . If  an  overlap  exists  between  any  two 
line  segments,  say  X^  and  X2»  then  X2  is  shortened  to  eliminate  the  gap. 

The  following  is  a description  of  the  inportant  variables  inRSS2. 

Label  Type  Description 

CONV  Real  Scaling  factor  to  convert  X-Y  points 

to  grid  resolution  units. 
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Label 

T^e 

Description 

DISK 

Integer 

Variable  which  defines  the  strip  feature 
output  device  unit.  Disk  = 9 

DX 

Real 

AX  for  any  consecutive  point  pair. 

DY 

Real 

A Y for  any  consecutive  point  pair 

IDDDX 

Integer 

Ax  value,  may  be  positive  or  negative 

IDDX 

Integer 

X displacement  of  line  segment  from  vertex. 

IDX p IDY 

Integer 

Integer  representations  of  DX,  DY 

IF 

Integer 

Feature  code  of  current  feature 

IFG 

Integer 

Flag  indicating  direction  of  line  seg- 
ment. Negative  value  means  left  edge, 
positive  value  means  right  edge,  zero 
means  horizontal. 

IFILL 

Integer 

Counter  flag  indicating  whether  fill 
strip  is  required. 

IGDSIZ 

Integer 

Grid  resolution  size  in  meters  X1000 

IIF 

Integer 

Merge  priority  code  for  feature 

I IT 

Integer 

Number  of  segments  in  the  JX+IX  arrays 

IN 

Integer 

Number  of  coordinates  in  feature  record 

INTAP 

Integer 

Variable  which  defines  the  input  device 
unit  number.  INTAP*=8 

IN2W 

Integer 

Variable  which  defines  the  output  device 
unit  number  for  2 word  buffer.  IN2W=20 

IRATN 

Integer 

Length  of  line  segment  required 

IREGSZ 

Integer 

Region  size  on  a side 

I REM 

Integer 

Same  as  TSEG 

IRES 

Integer 

Grid  points  per  line.  IRES=2048 

IRTND 

Integer 

X displacement  of  line  segment  from  vertex 

ISEG 

Integer 

Number  of  line  segments  in  feature  record 
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ISEQ 

Integer 

Sequence  number  of  feature 

ISET 

Integer 

Map  segment  counter 

ISPEC 

Integer 

Specularity  flag 

2 

IT 

Integer 

Hold  area  for  ISEG 

t 

ITH 

Integer 

Grid  scaled  theta  value 

i 

i 

ITRC 

Integer 

Counter  of  attempted  reads 

j 

ITYP 

Integer 

Feature  type.  l=point,  2=line, 
3= polygon 

i 

IX 

Integer 

Array 

Feature  output  buffer 

i 

i 

JSEG 

Integer 

Number  of  line  segments  in  JX  array 

! 

JT 

Integer 

Hold  area  for  JSEG 

JX 

Integer 

Array 

Array  of  line  segments  for  current 
strip 

j 

1 

KX 

Integer 

Array  of  grid  scaled  data  points 

i 

MWD 

Integer 

Number  of  words  read  for  this  feature 

\ 

NI 

Integer 

Top  of  IX  list 

~ * 

Nil 

Integer 

Number  of  entries  in  the  IX+JX  buffers 

! 

NIJ 

Integer 

Pointer  to  second  element  of  IX  list 

"i 

NJ 

Integer 

Top  of  JX  list 

J 

NK 

Integer 

Position  pointer  for  fill  strip  into 
IX  buffer 

i 

NSEG 

Integer 

Number  of  line  segments  required  to 
bridge  current  point  pair 

i 

NWD 

Integer 

Length  of  feature  record 

i 

PI 

Real 

Inverse  tangent  of  one  unit 

f 

RATIO 

Real 

Absolute  value  of  ^x/  AY 

t 

! 
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RATN 

Real 

Value  represents  X+(Ax/A^) 

RESK 

Real 

Resolution  element  size  in  feet 

SCALE 

Real 

Scale  factor.  (ex.  100000) 

TH 

Real 

Theta  value  unsealed  to  grid 

li 

WORD 

Real 

Array 

Array  for  2 word  output  buffer 

; 

X2 

Real 

Input  feature  buffer 

t .1 

Array 

X5 

Integer 

Table  for  conversion  of  feature  codes 

i 

Array 

merge  priority  codes 

2.2.2. 7 

Input. 

RSS2  requires 

one  disk  input  file  (Tape  8)  from  RSS1 

and  one 

card  input  file  with  the 

following  format. 

(i 

i 

Columns 

1-20 

Contains  the  conversion  factor  for 
changing  inches  on  the  map  to  feet 
on  the  ground.  For  exanfple,  if  the 
map  scale  is  1:1000,000,  this  field 
contains  the  number  8333.333 

li 

Columns 

21-35 

Blank 

l 

Columns 

36-40 

Contain  the  size  of  the  grid  resolution 
elements  used  to  label  points  in  a region 
of  planimetry  data.  This  value  must  be 
expressed  in  units  of  meters  X 1000. 

li 

Columns 

41-48 

Blank 

Oolumns 

49-50 

Contain  the  region  size,  i.e.  the  number 
of  resolution  elements  along  a region  edge. 
The  restrictions  on  this  value  are  (1)  it 
must  not  exceed  48  because  RSS7  which  pro- 

1 

i ; 

cesses  by  region  has  core  allocation  for  a 
maximum  38  x 48  region  size;  and  (2)  it 
must  be  a multiple  of  4 because  the  input 
data  for  RSS7  is  packed  four  data  items 
per  word  and  this  requirement  would  eliminate 

( * * 

any  need  for  unpacking  and  repacking  of  data. 

r 

u 

i 1 

j / i 
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2. 2. 2. 8 Output.  RSS2  outputs  two  disk  files.  The  first  output  file 
(Tape  9)  contains  1800  word  feature  records.  Each  record  has  a 6-word 
header  with  the  remaining  record  divided  into  4-word  strip  descriptions. 
The  format  of  the  feature  buffer  is  as  follows. 

Header  - Words  1-6 


Word 

1 

Feature  sequence  number 

Word 

2 

Feature  code 

Word 

3 

Number  of  coordinates  in  feature  buffer 

Word 

4 

Number  of  line  segments  in  feature  buffer 

Word 

5 

Specularity  flag  for  feature 

Word 

6 

Feature  type  l=point,  2=1 ine,  3=polygon 

Strip  Description  - Word  6-1800 

Word  6+4N-3  Y-value  for  the  Nth  strip 

Word  6+4N-2  X-value  for  the  beginning  of  the  Nth  strip 

Word  6-4N-1  AX-value  for  the  Nth  strip,  equals  the 

total  number  of  grid  points  covered  by 
the  strip.  The  endpoint  of  the  strip  is 
then  set  to  X+  X-l. 


Word  6+4N 


Specularity  angle  for  the  Nth  strip.  This 
quantity  is  used  to  determine  abnormal 
radar  return  quantities.  Presently  not  used. 


The  second  output  file  (Tape  20)  consists  of  a single  two-word  entry. 

Word  1 Resolution  element  size  in  meters  X 1000 

Word  2 Number  of  resolution  elements  along  a 

region  edge 

This  file  is  created  on  disk  to  eliminate  the  need  for  card  input  of  the 
same  data  throughout  the  remainder  of  the  system. 

2. 2. 2. 9 Externals.  System  routines  called  by  RSS2  are  LENGTH  and  UNIT 
which  are  mass  storage  1/0  function  routines,  and  ATAN2  which  computes  the 
inverse  tangent  between  any  two  arguments.  There  are  no  subroutine  calls 
in  this  program. 
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2.2.2.10  Error  Conditions.  RSS2  has  three  error  conditions.  One  error 
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condition  occurs  if  the  card  input  resolution  size  or  region  size  is  out 
of  the  allowable  range.  The  error  messages  output  for  this  condition  are 
as  follows: 

INPUT  ERROR- RESOLUTION  ELEMENT  TOO  FINE 
OR 

INPUT  ERROR-REGION  SIZE  TOO  LARGE 
Another  error  condition  occurs  if  there  is  an  error  found  while 
writing  either  of  the  two  output  buffers  of  this  program.  The  error 
message  output  is  as  follows: 

OUTPUT  ERROR  ON  UNIT  nn 

The  next  error  condition  occurs  if  an  error  is  found  while  read- 
ing the  feature  input  buffer.  The  message  for  this  condition  is  as  follows: 

IN  FEATURE  NO.  nnnn  ERROR  AT  Y=nnnn 

2.2.2.11  Program  Flowchart.  The  flow  diagram  for  RSS2  is  illustrated  in 
Figure  I I. 4 

2.2.2.12  Program  Listing.  The  program  listing  for  RSS2  is  attached. 

2.2.3  and  2.2.4 

2. 2. 3.1  Program  Name.  RSS3 

2. 2. 3. 2 Storage . Approximately  21K  words  of  core 

2. 2. 3. 3 Run  Time.  Approximately  230  cp  seconds  and  2 minutes  wall 
clock  time. 

2. 2. 3.4  Cost.  Approximately  35  dollars 

2. 2. 3. 5 Program  Function.  RSS3  assigns  the  planimetry  data  strips  to 
regions.  If  a strip  overlaps  several  regions,  say  N of  them,  it  is  broken 
up  into  N segments  with  each  segment  being  assigned  to  its  respective  region. 
RSS3  also  assigns  merge  priorities  for  each  feature  code. 

2. 2. 3.6  Program  Description.  The  procedure  used  to  assign  the  planimetry 
strips  to  their  respective  regions  is  a simple  one.  Consider  for  example 
a strip  with  the  following  characteristics: 

2-25 


'4b  . * . ■ ...  x ' 1 1 1 iltiw— OlWMisui  I 

__ 


i 


t 


Y = 750 
X = 875 
Ax  = 24 

Region  Size  = 32 
Resolution  = 156.25  feet 

This  means  that  the  strip  begins  117,187  feet  north  (750  X 156.25) 
and  136,718  feet  east  (875  X 156.25)  of  the  map  origin.  We  have  assumed 
that  each  region  contains  32  increments  along  each  edge  so: 

*Y  = 750/32  = 23.43750 

R = 875/32  = 27.34375 
x 

If  we  assume  that  each  horizontal  row  of  regions  can  contain  at  most  331 
regions,  then  the  region  number  for  the  beginning  of  the  strip  is: 

Nr  = 331  * (23)  + 28  = 7641 

The  (x, Y)  coordinates  of  the  beginning  point  relative  to  the  first  grid 
element  in  the  region  is: 

Y * 0.43750  * (32)  = 14 
X = 0.34375  * (32)  = 11 

Since  X = 24,  the  end  point  of  the  strip  is  at  X = 11+24-1  = 34. 

The  maximum  address  within  a region  is  32,  so  this  strip  overflows  into 


the  next  region 

to  the  east.  We 

therefore  wind 

up  with 

two  strips  as 

follows: 

STRIP  1 

STRIP  2 

Region 

7641 

Region 

7642 

Y 

14 

14 

X 

11 

1 

A X 

22 

2 
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Numbering  of  these  regions  begins  in  the  southwest  corner  of  the 
map  with  the  assignments  being  made  sequentially  left  to  right  and  bottom 
to  top.  The  merge  priorities  assigned  to  each  feature  are  done  through  the 
use  of  a translation  table  which  is  hard-coded  as  a DATA  statement. 

The  following  is  a description  of  the  important  variables  in  RSS3. 


Label 

Type 

Description 

DISK 

Integer 

Variable  which  defines  the  input  unit 
device  number  for  the  planimetry  strip 
feature  data.  DISK-7 

IF 

Integer 

Feature  code 

IFC 

Integer 

Array 

Merge  priority  conversion  table 

II 

Integer 

Y value  for  the  southwest  corner  of 
the  region 

III 

Integer 

Region  count  in  Y direction 

ILIM 

Integer 

Pointer  to  segment  count  in  feature 
record. 

I MAX 

Integer 

Maximum  X value  for  Y 

IMIN 

Integer 

Minimum  X value  for  Y 

IN 

Integer 

Array 

1800  word  input  buffer  for  strip 
feature  record 

INRC 

Integer 

Number  of  records  output 

INSG 

Integer 

Number  of  segments  output 

IN2W 

Integer 

Variable  which  defines  the  input  device 
number  for  the  2-word  file  from  RSS2. 
IN2W=20 

I REG 

Integer 

Region  number 

IRES 

Integer 

Number  of  points  per  plot  line 

ISEG 

Integer 

Number  of  segments  in  record 
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ISEQ 

Integer 

Sequence  number 

ISPEC 

Integer 

Specular  code 

1ST 

Integer 

Counter  for  number  of  words  of  feature 
data  input 

ITYP 

Integer 

Feature  type 

I XL1 

Integer 

Region  number  for  IMIN 

IXL2 

Integer 

Region  number  for  IMAX 

IYL1 

Integer 

Y value  for  bottom  of  feature 

IYL2 

Integer 

Y value  for  top  of  feature 

JDX 

Integer 

Array 

Array  of  Ax  values  for  disk  output 

JJ 

Integer 

X value  for  the  southwest  corner  of 
the  region 

JJJ 

Integer 

Region  count  in  X direction 

JTH 

Integer 

Array 

Array  for  8 values  for  disk  output 

JX 

Integer 

Array 

Array  for  X values  for  disk  output 

JY 

Integer 

Array 

Array  for  Y values  for  disk  output 

MWD 

Integer 

Number  of  words  transferred  from 
input  buffer 

NRX 

Integer 

Number  of  regions  in  the  X direction 

NVM 

Integer 

Number  of  segments  in  record 

NX 

Integer 

Number  of  elements  in  X direction  of 
region 

NY 

Integer 

Number  of  elements  in  Y direction  of 
region 

OUTAP 

Integer 

Variable  which  defines  output  unit 
device  number.  OUTAP=l 

WORD 

Integer 

Array 

Array  for  2-word  buffer  output  by  RLMS2 
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2. 2. 3. 7 Input.  RSS3  requires  two  disk  input  files.  The  files  required 
are  the  two  word  output  file  (Tape  20)  from  RSS2,  and  the  strip  feature 
data  file  (Tape  7)  from  RSS2. 

2. 2. 3.8  Output.  RSS3  outputs  one  disk  file  (Tape  2)  in  the  form  of  card 
images.  The  format  for  each  image  is  as  follows: 


Word  1-5 
Word  6-7 
Word  8-9 
Word  10 
Word  11 
Word  12 
Word  13-89 


Region  number 

Merge  priority  code 

Feature  code 

Feature  type 

Specularity  code 

Number  of  segments  in  record 

Y , X, Ax,  0 for  up  to  seven  segments 


The  only  printout  from  RSS3  is  a message  indicating  the  end  of  processing. 


2. 2. 3.9  Externals.  System  routines  called  by  RSS3  are  LENGTH  and  UNIT 
which  are  mass  storage  I/O  function  routines.  There  are  no  subroutine 
calls  in  this  program. 

2.2.3.10  Error  Conditions.  RSS3  has  two  error  conditions.  One  error 
condition  occurs  when  there  is  an  error  found  while  reading  the  two  word 
disk  file  output  by  RSS2.  The  error  message  for  this  condition  is  as 
follows: 


INPUT  ERROR  ON  UNIT  20 

The  other  error  condition  occurs  when  an  error  is  found  while  reading  the 
input  strip  feature  file  output  from  RSS2.  The  error  message  for  this 
condition  is  as  follows: 

BUFFER  IN  ERROR  AT  ISEQ  = XXXXXXX 

2.2.3.11  Program  Flowchart.  The  flow  diagram  for  RSS3  is  illustrated 
in  Figure  II.  5 

2.2.3.12  Program  Listing.  The  program  listing  for  RSS3  is  attached. 


I 

I 

I 
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I Read  File 
' Output  j 
From  RLMS2/ 


X-RMINX 


Read 

Feature 

Record 


Sort  Re- 
cords By 
Region, 
Merge  Pri- 


J=J+1 


fine  Region 
Sounds  For  Top 
Snd  Bottom  of 
Feature 


Begin  ir 
^Region  xy 
vRpw 


Y=YMAX 


>/^tri?V 

Extend  X 
^cross  Re^ 
^gion  Bpds 


Assign 
Strip  to 
Region 


i 


Y < YMIN 


Divide  Strip! 
at  Region  | 
Boundary  ,1 


Determine 
MINX  and 
MAXX  for 
Y 


Assign 
Section  on 
Left  to 
Region 


Calculate 
Region  Bounds 
( RMINX , RMINY) 
?or  MINX  and 


Define  New 
Strip  From 
Right  Hand 
Section 


FIGURE  II. 5 - PROGRAM  RSS3  FLOWCHART 
(Page  1 of  2)  “ 
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2. 2.4.1  Program  Name.  RSS4  (Sort  Routine) 

2. 2.4.2  Program  Function.  RSS4  is  a sorting  routine  which  reorders 
the  planimetry  data  generated  by  RSS3  by  (1)  ascending  region  number; 

(2)  ascending  merge  priority;  (3)  increasing  Y value;  and  (4)  increasing 
X value. 

2. 2. 4. 3 Program  Description.  RSS3  employs  the  CDC  SORT/MERGE  package 
to  order  the  card  images.  For  the  purposes  of  this  document,  we  will 
call  this  routine  RSS4. 

In  preparation  for  RSS5,  the  card  images  are  sorted  in  the 
following  order; 


ORDER  ATTRIBUTE 

1 Region  Number 

2 Merge  Priority 

3 Y-value  of  first  strip 

4 X-value  of  first  strip 

Although  the  reader  is  directed  to  the  CDC  SORT/MERGE  manual  for 

conqplete  details  on  the  operation  of  the  SORT/MERGE  package,  we  present 

here  a listing  and  brief  description  of  RSS4. 

SORT,  VAR=DISK 
BYTESIZE,6 

FILE, SORT=TAPE 1 , OUTPUT=TAPE  2 

FIELD , REGION (1.1,5, DISPLAY) , PRIORITY (6.1,2, DISPLAY) , 

, INTY  (13. 1,2, DISPLAY)  , INTX (15. 1 , 2 , DISPLAY) 

KEY, REG ION (A, OWN)  .PRIORITY  (A, OWN) , INTY (A, OWN) , INTX(Z .OWN) 
SEQUENCE , OWN ()rf, 0,1, 2, 3, 4, 5, 6, 7, 8, 9) 

OPTIONS , RETAIN 
END 

The  FIELD  card  lists  the  attributes  to  be  sorted  in  order  of 
decreasing  priority.  It  indicates  the  format  of  the  data,  which  in  our 
case  is  CDC  DISPLAY  code,  and  its  location  within  the  record.  The  KEY 
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card  gives  the  direction  of  the  sort  (i.e.  the  "A"  indicates  that  the 
numbers  are  to  be  placed  in  ascending  order)  and  the  ordering  scheme,  which 
in  our  case  is  the  sequence  OWN.  The  SEQUENCE  card  is  used  to  define  an 
ordering  sequence  other  them  one  of  the  standard  sequences.  This  card  is 
required  in  RSS4  because  a blank  in  display  code  is  represented  by  55B 
and  this  value  is  greater  than  that  assigned  to  the  integers.  Without  a 
special  sequence,  a number  like  W1023  would,  for  example,  be  placed  after 
1416486  in  an  ascending  order  sort  of  the  card  images. 

Clearly,  the  output  from  ASS 4 is  in  the  same  format  as  that  from 
RSS3.  Only  the  ordering  of  the  records  is  changed. 

2. 2.4.4  Input.  RSS4  requires  the  disk  output  planimetry  data  file 
(Tape  2)  from  RSS3. 

2.2.4. 5 Output . The  disk  output  from  RLMS4  is  the  same  format  as  that 
from  RSS3.  Only  the  ordering  of  the  records  has  changed. 

2.2.5 

2. 2. 5.1  Program  Name  RSS5 

2. 2. 5. 2 Storage  Approximately  44K  words  of  core 

2. 2.5. 3 Run  Time  Approximately  217  cp  seconds  and  2 minutes  wall 
clock  time 

2. 2. 5. 4 Cost.  Approximately  45  dollars 

2. 2. 5. 5 Program  Function.  RSS5  merges  the  planimetry  strips  so  that  no  strips 
overlap.  It  marks  the  end  of  the  planimetry  data  preparation. 

2. 2. 5.6  Program  Description.  The  records  containing  the  planimetry  strips 
have  already  been  sorted  by  region  number  and  merge  priority.  The  merge 
priority  scheme  is  set  up  so  that  features  which  may  be  overwritten  appear 
before  those  which  may  not.  For  example,  consider  two  features,  a city 
(low  merge  priority)  and  a lake  (higher  merge  priority) . Each  of  these  has 
associated  with  it  a series  of  strips  describing  its  shape  and  location. 
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Now,  if  the  lake  lies  within  the  bounds  of  the  city,  a merge  of  the  data 
is  required  since  it  is  necessary  to  specify  which  feature  is  to  be  "on 
top";  i.e.  the  lake  strips  must  be  copied  over  the  city  strips  or  the  lake 
will  not  appear  on  the  final  picture.  The  scheme  is  presently  imperfect 
(consider  an  island  with  a lake  on  it)  but  can  be  expected  to  work  in  the 
vast  majority  of  cases. 

In  essence,  rss5  simply  copies  the  planimetry  strips  into  a core 
array  in  the  order  in  which  they  appear  on  the  input  file.  This  array 
becomes  an  image  of  the  region  structure  as  successive  strips  are  copied 
into  it.  An  enpty  record  is  written  for  those  regions  without  planimetry 
data.  As  an  example  of  the  operation  of  RSS5,  consider  two  planimetry 
strips,  one  being  part  of  a lake  and  the  other  being  part  of  a city.  Take 
the  location  and  size  of  these  strips  to  be  as  follows: 

Lake  Strip  City  Strip 


Y = 25 
X = 11 

Ax  = 6 


Y = 25 
X = 8 

A*  = 15 


Clearly,  this  portion  of  the  lake  overlaps  the  city.  RSS5  takes  these 
strips  and  generates  three  (3)  strips,  two  belonging  to  the  city  and  one 
belonging  to  the  lake.  The  location  and  size  of  these  strips  is  as  follows: 
Lake  Strip  City  Strip 


Y = 25 
X = 11 

A x * 6 


Y = 25 
X * 8 

Ax  = 3 


Y = 25 
X = 17 

Ax  = 6 


Notice  that  these  strips  do  not  overlap. 
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The  following  is  a description  of  the  important  variables  in  RSS5 


Label 

DISK 

IEND 

IF 

IG 

IKK 

IKC 

IKC2 

IN 

INDEX 

IN2W 

I REG 
I SUM 

ITH 

IX 

IY 

IZ 


Ty£g 

Integer 


Integer 

Integer 

Array 

Integer 

Array 

Integer 

Integer 

Integer 

Integer 

Integer 

Array 

Integer 


Integer 

Integer 


Integer 

Array 

Integer 

Array 

Integer 

Array 

Integer 

Array 


Description 

Variable  which  defines  the  input  unit 
device  number  for  the  planimetry  data. 
DISK=7 

End  of  file  flag 

Common  block  of  feature  codes  for  first 
or  previous  strip 

Common  block  of  merge  priority  codes 
for  first  or  previous  strip 

Merge  priority  for  input  strip 

Region  size  of  planimetry  data  on  a side 

Total  number  of  elements  in  a region 

Number  of  segments  input 

Storage  area  required  by  mass  storage 
I/O  routines 

Variable  which  defines  the  input  device 
number  for  the  2 word  file  from  RSS2 
IN2W=20 

Input  region  number 

Total  number  of  elements  in  the  current 
region 

Common  block  of  specularity  angles  for 
first  or  previous  strip 

Common  block  of  Y-values  for  first  or 
previous  strip 

Common  block  of  X-values  for  first  or 
previous  strip 

Common  block  of  Ax  values  for  first  or 
previous  strip 
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JF 

Integer 

Array 

Common  block  of  feature  codes  for 
current  strip 

JG 

Integer 

Array 

Conraon  block  of  merge  priority  codes 
for  current  strip 

JN 

Integer 

Number  of  segments  in  current  strip 

JREG 

Integer 

Current  region  number 

JTH 

Integer 

Array 

Common  block  of  specularity  angles 
for  current  strip 

JX 

Integer 

Array 

Common  block  of  Y values  for  current 
strip 

JY 

Integer 

Array 

Common  block  of  X values  for  current 
strip 

JZ 

Integer 

Array 

Common  block  of  Ax  values  for  current 
strip 

KF 

Integer 

Array 

Common  block  work  area  for  merged  feature 
codes 

KG 

Integer 

Array 

Common  block  work  area  for  merged  priority 
codes 

KTO 

Integer 

Array 

Common  block  work  area  for  merged  specu- 
larity angles 

KX 

Integer 

Array 

Compton  block  work  area  for  merged  Y values 

KY 

Integer 

Array 

Common  block  work  area  for  merged  X values 

KZ 

Integer 

Array 

Common  block  work  area  for  merged  &X 
values 

N,K 

Integer 

Number  of  elements  (words)  in  region 

NBYT 

Integer 

Counter  for  number  of  elements  output 

NEMT 

Integer 

Counter  for  number  of  empty  regions  output 

NFUL 

Integer 

Counter  for  number  of  filled  regions  output 

NSEG 

Integer 

Counter,  not  used 
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NTOT 

Integer 

Counter  not  used 

NY 

Integer 

Maximum  number  of  regions 
direction 

in  Y 

NX 

Integer 

Maximum  number  of  regions 
direction 

in  X 

WORD 

Integer 

Array  for  two  word  buffer  output  by 
RLMS2 

2. 2. 5. 7.  Input  RSS5  requires  two  input  disc  files.  The  files  required 
are  the  two  word  output  file  (TAPE  20)  from  RSS2,  and  the  sorted  output 
file  (TAPE  7)  from  RSS4. 

2. 2. 5. 8 Output.  RSS5  outputs  one  random  access  disc  file  (TAPE  10). 

This  output  file  consists  of  a record  for  each  region  on  the  map 
regardless  of  whether  or  not  that  region  actually  contains  any  planimetry 
data.  The  first  word  of  the  record  contains  the  total  number  of  strips 
contained  in  the  region,  while  the  remaining  words  contain  the  information 
describing  the  strips.  Each  strip  is  described  by  one  word,  with  the 
data  being  packed  as  follows: 


BITS 

31  - 36 
25  - 30 
19  - 24 
7-18 
1-6 


The  only  printout  from  RSS5  is  a 


INFORMATION 

Y - coordinate  of  strip 
X - coordinate  of  strip 
AX  for  strip 

Specularity  angle  for  scrip 
Feature  code  for  strip 

isage  indicating  the  end  of  processing. 


2. 2. 5.9  Externals . System  routines  called  by  RSS5  jure  OPENMS,  CLOSMS, 
WRITMS , and  UNIT  which  are  mass  storage  I/O  function  routines.  The 
subroutine  called  by  RSS5  is  RADOUT.  This  subroutine  is  discussed  under 
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section  2.2.5.13 


2.2.5.10  Error  Conditions.  RSS5  has  three  error  conditions.  One 


error  condition  occurs  if  there  is  an  error  found  while  reading  the 
two  word  disc  file  output  by  RSS2.  The  error  message  for  this  condition 
is  as  follows: 

INPUT  ERROR  ON  UNIT  20 

Another  error  condition  occurs  if  the  total  number  of  elements  ISUM 
output  by  a region  is  not  equal  to  the  square  of  the  region  size  on  a 
side  (IKC2).  If  ISUM  is  larger  than  IKC2,  then  the  500  word  common 
blocks  will  overflow.  If  ISUM  is  smaller  than  IKC2,  part  of  the  region 
record  will  be  empty.  The  error  message  output  for  this  condition  is 
as  follows: 

PARITY  ERROR  DURING  INDEX  WRITE-OUT 

2.2.5.11  Program  Flowchart.  The  flow  diagram  for  RSS5  is  illustrated  in 
Figure  II. 6. 

2.2.5.12  Program  Listing.  The  program  listing  for  RSS5  is  attached. 

2.2.5.13  Subroutine  Description. 

2.2.5.13.1  Subroutine  Name.  RADOUT 

2.2.5.13.2  Summary ■ RADOUT  packs  and  writes  out  the  planimetry  region 
records . 

2.2.5.13.3  Description  of  Processing.  The  calling  sequence  for  RADOUT 
is  - RADOUT (K) . K is  the  number  of  features  in  the  region.  For  each 
feature,  pertinent  information  is  packed  in  to  one  word  describing  that 
feature  and  stored  into  a region  record.  When  the  region  record  packing 
is  complete,  the  record  is  output  to  a random  access  file.  For  a 


2-40 


FIGURE  II. 6 - PROGRAM  RSS5  FLOWCHART 
(Page  1 of  4) 
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description  of  the  important  variables  in  this  subroutine  refer  to  the 


list  under  section  2. 2. 5.6. 

2.2.5.13.4  Error  Conditions.  None 

2.2.5.13.5  Subroutine  Flowchart.  A flow  diagram  is  not  considered  necessary _ 
Consult  the  listing  of  RSS5  for  details. 

2.2.5.13.6  Subroutine  Listing.  The  listing  for  RADOUT  is  included  with 
that  of  RSS5. 

2.2.6 

2. 2. 6.1  Program  Name.  RSS6 

2. 2. 6. 2 Storage.  Approximately  64K  words  of  core. 

2.2.6. 3 Run  Time.  Approximately  117  cp  seconds  and  1 minute 
wall  clock  time. 

2. 2. 6. 4 Cost.  Approximately  43  dollars. 

2. 2. 6. 5 Program  Function.  RSS6  divides  the  terrain  data  into  region 
blocks . 

2. 2. 6. 6 Program  Description.  The  input  data  base  is  characterized  by 
resolution  elements  whose  length  is  defined  by  word  34  of  the  header  record, 
in  units  of  meters  x 1000.  The  resolution  elements  are  arranged  in  horizontal 
strips  covering  the  entire  width  (east-west  dimension)  of  the  map  sheet. 

Word  two  of  a two  word  disc  input  file  defines  the  region  size  in  resolution 
elements  or  a side.  Horizontal  strips  of  data  whose  length  is  defined  by 
the  region  size  is  read  into  core.  This  data  is  reformatted  into  rows  of 
region  blocks  numbered  from  the  lower  left  corner  of  the  map.  Region 
numbers  are  determined  by  words  21  and  22  of  the  header  record.  Word  21 
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defines  the  number  of  regions  in  the  Y direction  and  word  22  defines  the 
number  of  regions  in  the  X direction.  After  all  the  regions  of  the  map 
formatted  and  output  to  a temporary  random  disk  file,  the  region  records 
are  reread  in  order  and  output  to  a permanent  sequential  disk  file. 

The  following  is  a description  of  the  important  variable  for  RLMS6. 


LABEL 

TYPE 

DESCRIPTION 

ADDR 

Integer 

Pointer  to  starting  address  of  data 
on  mass  storage  device. 

I ARRAY 

Integer 

Array 

Temporary  holding  area  of  data  for  one 
region. 

IK 

Integer 

Region  size  of  data  (in  meters  x 1000) 

IM 

Integer 

Number  of  words  in  a region  record. 

IN 

Integer 

Number  of  words  required  to  read  in 
data  for  one  row  of  regions. 

INBUF 

Integer 

Array 

Input  buffer  for  terrain  tape  data. 

INDEX 

Integer 

Array 

Storage  area  required  by  the  mass 
storage  I/O  routines. 

INTAP 

Integer 

Variable  which  defines  the  input 
unit  device  number  for  the  terrain 
data.  INTAP  * 8 

INZW 

Integer 

Variable  which  defines  the  input 
device  number  for  the  2 word  file  from 
RSS2.  INZW  = 20. 

I POINT 

Integer 

Pointer  to  data  item  within  current  row 
of  regions. 

JJ 

Integer 

Region  row  number 

K 

Integer 

Number  of  words  currently  in  region 
record  currently  being  formatted. 

MOVE 15 

Integer 

Array 

Array  defining  the  number  of  bytes  to 
shift  a word  of  the  packed  input  record 

NCOL 

Integer 

Number  of  data  elements  in  the  X 
direction. 

NREC 

Integer 

Number  of  records 

N REGION 

Integer 

Region  number 

NROW 

Integer 

Number  of  data  elements  in  the  Y 
direction. 

NXREG 

Integer 

Number  of  regions  in  the  X direction. 

NY  REG 

Integer 

Number  of  regions  in  the  Y direction. 

OUT 

Integer 

Array 

Buffer  for  formatted  output  region 
record. 

OUTAP 

Integer 

Variable  which  defines  output  device 
unit  number  for  formatted  region  record 
OUTAP  - 9 

WORD 

Integer 

Array 

Array  for  two  word  buffer  output  by 
RSS2. 

XCOOR 

Integer 

Current  X coordinate  position. 

XFILL 

Integer 

Number  of  filler  points  required,  if 
any,  in  the  X direction. 

YCOOR 

Integer 

Current  Y coordinate  position. 

YFILL 

Integer 

Number  of  filler  points  required,  if 
any,  in  the  Y direction. 

2 . 2 . 6 . 7 Input. 

RSS6  requires 

one  input  tape  file  and  one  input  disk 

file.  The  tape 

file  (TAPE  8)  contains  the  terrain  data  base.  The  disk 

file  (TAPE  20) 

contains  the  two 

word  buffer  output  by  RSS2. 

The  first  record  on  the  tape  is  a 36  word  header.  Each  entry  in 
the  header  is  an  integer  representation  of  the  actual  valve  multiplied 
by  1000  and  truncated.  The  following  is  the  description  of  the  header. 
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Word  21:  The  total  number  of  records  (profiles)  on  the  tape, 

not  counting  the  header. 

Word  22:  The  total  number  of  elevation  points  per  record. 

Elevations  are  unsigned  15  bit  integers,  in  feet, 
packed  4 per  word. 

Word  23-25: 

Geographic  latitude  in  degrees,  minutes,  and  seconds 
of  the  radar  target. 

Word  26:  Ground  distance  in  meters  between  the  UNAMACE  point 

of  tangency  and  the  radar  target. 
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Word  27-29: 

Hie  longitude  in  degrees,  minutes,  and  seconds  of  the 
radar  target. 

Word  30:  The  ground  distance  in  meters  from  the  UNAMACE  point 

of  tangency  west  to  the  target. 

Word  31:  Hie  ground  distance  in  meters  from  the  first  elevation 

profile  (upper  NW  corner  of  the  elevation  data)  north 
to  the  radar  target.  A negative  value  indicates  that 
the  target  is  south  of  the  first  elevation  point. 

Word  32:  Hie  ground  distance  in  meters  from  the  first  elevation 

profile  west  to  the  target.  Hiis  number  is  also  negative. 

Word  33:  Target  ID  Code 

Word  34:  Ground  spacing  between  points  of  a profile  in  meters. 

Presently  this  value  must  be  equal  to  the  grid  resolu- 
tion size  output  by  RSS2. 

Word  35:  Ground  spacing  between  profiles  in  meters.  Equals 

word  34  at  present  time. 

Word  36:  Zero  fill  word. 


Following  the  header  are  a series  of  records  - one  for  each  profile 
line.  Each  profile  contains  the  elevation  values  for  a west-to-east  run- 
ning strip  of  terrain,  with  the  first  profile  corresponding  to  the  northern- 
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2. 2.6. 8 Output  RSS6  has  one  disk  output  file  (Tape  9)  arranged  such 
that  the  information  for  each  region  is  stored  in  a one-dimensional  array. 

The  first  word  contains  the  region  number  while  the  remaining  words  contain 
the  elevation  data  packed  four  per  word.  The  data  is  arranged  such  that 
the  southwestern  edge  of  the  map  appears  first.  The  only  printout  from 
RSS6  is  a message  indicating  the  end  of  processing. 

2. 2.6.9  Externals  System  routines  called  by  RSS6  are  OPENMS,  READMS , 
WRITMS,  and  UNIT  which  are  mass  storage  I/O  function  routines.  There  sure 
no  subroutine  calls  made  by  this  program. 

2.2.6.10  Error  Conditions.  RSS6  has  two  error  conditions.  The  first 
error  condition  occurs  if  there  is  an  error  during  input  or  output  to  the 
disk  files.  One  of  the  following  messages  is  output  for  this  error  condition 

PREMATURE  EOF  ENCOUNTERED  ON  INPUT  TAPE 
OR 

PARITY  ERROR  ENCOUNTERED ERROR  EXIT 

The  second  error  condition  occurs  if  the  resolution  value  for  the  terrain 
data  is  different  from  that  for  the  planimetry  data.  The  message  output 
for  this  condition  is  as  follows: 

ERROR-CULTURE  RESOLUTION  NOT  EQUAL  TO  TERRAIN  RESOLUTION 

2.2.6.11  Program  Flowchart.  The  flow  diagram  for  RSS6  is  illustrated  in 
Figure  II.  7. 

2.2.6.12  Program  Listing.  The  program  listing  for  RSS6  is  attached. 

2.2.7 

2. 2. 7.1  Program  Name.  RSS7 

2. 2. 7. 2 Storage . Approximately  42K  words  of  core 
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2.2.7. 3 Run  Time.  Approximately  36  cp  seconds  and  2 minutes  wall 
clock  time. 

2. 2. 7. 4 Cost.  Approximately  11  dollars. 

2. 2.7. 5 Program  Function.  Program  RSS7  merges  the  planimetry  file  from 
RSS5  and  terrain  file  from  RSS6  for  all  regions  lying  within  the  radar 
ground  range  of  the  target. 

2. 2. 7. 6 Program  Description.  Using  information  from  the  terrain  file 
header  record  and  the  card  input  record,  the  region  number  and  X-Y  coord- 
inates of  the  radar  target  are  determined.  The  target  location  coordinate 
values  are  relative  to  the  southwest  corner.  With  the  target  point  as 
its  center,  an  imaginary  circle  is  then  constructed  whose  radius  is  equal 
the  ground  range  of  the  radar.  A sqaure  circumscribed  about  this  circle 
will  then  contain  all  of  the  regions  required  to  produce  the  radar  scene 
at  the  chosen  altitude;  the  rest  may  be  discarded.  This  procedure  is 
illustrated  in  Figure  II.  10.  For  each  region  overlapped  by  the  circum- 
scribed square,  beginning  at  the  southwest  corner,  the  terrain  and  plani- 
metry data  are  merged  into  a single  record  labeled  by  the  region  number  and 
output  to  a disk  file.  Before  termination  a parameter  file  with  related 
map  quantities  needed  to  construct  the  radar  scene  is  set  up  and  output 

to  a disk  file.  The  list  of  important  variables  for  this  program  is  as 
follows : 

Label  Type  Description 

ALTTDE  Real  The  height  in  feet  of  the  radar  above  the 

target.  Maximum  value  allowed  is  32000. 
Default  value  is  32000.  This  variable  is 
part  of  the  parameter  list. 


I 


I 


1 

Label 

Type 

Description 

DEPNGL 

Real 

The  radar  depression  angle  function  which 
is  equal  to  1.5  + cot  25°  = 3.6445. 

DATA 

Integer 

Array 

Array  for  500  word  input  culture  record  and 
variable  length  output  merged  record. 

1 

DLINK 

Real 

Array 

Array  for  parameter  output  record  of  Real 
data  values. 

1 

ICULT 

Integer 

Switch  for  culture  data.  1=  culture  data, 
o=no  culture  data.  Default  is  one. 

■ 

1 

IK 

Integer 

Region  size.  This  variable  is  set  equal  to 
the  value  input  in  RSS2.  Ibis  variable 
must  be  a multiple  of  4 and  have  a maximum 
value  of  48.  This  is  part  of  the  parameter 
list. 

1 

\ m 

IFREQ 

Integer 

The  number  of  radial  scan  lines  to  be  con- 
structed per  degree.  Ibis  variable  is  part 
of  the  parameter  list. 

1 

■ 

ILIM,  JLIM 

Integer 

Computer  column  and  row  of  upper  corner 
region  of  radar  range. 

1 

IM 

Integer 

Number  of  words  in  terrain  data  input  record 

1 

INDEX, 

IDEXM1 

Integer 

Pointers  into  LDATA  array 

1 

IN2W 

Integer 

Input  buffer  for  two  word  file  output  by 
RSS2. 

| 

IPNT 

Integer 

Pointer  used  to  isolate  the  target  altitude 
value  in  the  packed  output  buffer. 

JB 

I REG, REG 

Integer 

Current  region  number  being  processed. 

i 9 

ISHIPT 

Integer 

Computed  number  of  bits  to  shift  word. 

I i 

ITALT 

Integer 

Target  altitude  in  feet  above  sea  level. 

i 

I TERR 

Integer 

Switch  for  terrain  data.  l>terrain  data. 
o=no  terrain  data.  Default  is  one.  Ibis 
variable  is  part  of  the  parameter  list. 

ITREG 

Integer 

Region  number  of  target 

I 


Label 

Type 

Description 

ITX.ITY 

Integer 

Computed  number  of  regions  in  X and  Y 
directions 

IX,  IY 

Integer 

Computed  column  and  row  of  lower  corner 
region  of  radar  range 

IXCOR, 

IYCOR 

Integer 

Location  of  target  within  region 

Kl,  K2 

Integer 

Pointers  into  the  culture  data  buffer 

LDATA 

Integer 

Array 

Array  for  parameter  output  record  of 
integer  data  values 

MOVE 

Integer 

Array 

Number  of  bits  to  shift  word 

MSIX 

Integer 

Array 

Index  directories  used  by  mass  storage  I/O 

NBASER 

Integer 

Starting  region  number  of  map  relative  to 
the  southwest  corner 

NR 

Integer 

Number  of  regions  of  radar  coverage 

NX 

Integer 

Maximum  number  of  regions  in  X direction. 
This  variable  is  part  of  the  parameter  list. 

RANGE 

Real 

Radar  radius  of  coverage  in  nautical  miles. 
This  variable  is  part  of  the  parameter  list 

RES 

Real 

Resolution  unit  size  (in  feet) 

REGNM 

Real 

Regions  per  nautical  mile 

REST 

Real 

Resolution  unit  size  (in  meters  X 1000) . 
This  variable  is  part  of  the  parameter  list 

TARGALT 

Real 

Target  altitude  in  feet.  This  variable  is 
part  of  the  parameter  list 

TABLE 

Real 

Array 

Array  for  input  terrain  record  packed  four 
heights  per  word. 

X,  Y 

Real 

Radar  location  in  meters  from  SH  origin. 
These  variables  are  part  of  the  parameter 
list. 

Label 


T*E£ 


Description 


Integer 

Array 


Array  for  two  word  file  output  from  RSS2. 


XLON, 

YLAT 


X-Y  target  coordinates  computed  relative  to 
the  northwest  corner. 


YDIST 


Distance  in  Y-direction  of  the  target  in 
nautical  miles  from  the  northwest  corner 
of  the  map. 


2. 2. 7.  7 Input.  RSS7  requires  three  disk  files  and  one  card  input  file. 

The  disk  files  required  are  the  output  terrain  file  (Tape  1)  from  RSS6,  the 
output  culture  file  (Tape  3)  from  RSS5,  and  the  two  word  file  (Tape  20) 
output  from  RSS2.  The  card  input  supplies  certain  NAMELIST  quantities.  The 
NAMELIST  variables  are  as  follows: 

I FREQ,  I TERR,  ALTTDE,  X,  Y 

See  the  program  description  for  details  of  these  variables. 

2. 2. 7. 8 Output . RSS7  outputs  two  disk  files.  The  first  file  (Tape  6)  is 
a parameter  file.  It  contains  two  five-word  records  containing  the  descrip- 
tive information  required  in  the  construciton  of  the  radar  scene.  Record 
one  contains  the  following  REAL  variable  data: 


Word  1 
Word  2 
Word  3 
Word  4 
Word  5 


RANGE 

ALTTDE 

X 

Y 

TARGALT 


Record  two  contains  the  following  INTEGER  variable  data: 


Word  1 
Word  2 
Word  3 
Word  4 
Word  5 


NX 

REST 

IFREQ 

IK 

ITERR 


See  the  program  description  for  details  of  these  variables. 
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The  second  file  (Tape  12)  is  a randoin-access  file  containing  merged  plani- 
metry-terrain  data  records.  These  records  are  labeled  by  region  number  and 
have  the  following  format: 

Word  1 Region  number 

Word  2 Descriptive  information  for  the  N strips 

contained  in  the  region.  Each  word  con- 
tains the  information  for  one  strip,  as 
output  by  RSS5. 

Packed  elevation  data  as  output  by  RSS7. 

A total  of  (IK  /4+1)  words  are  required 
to  store  the  (IK  X IK)  elevation  values 
for  each  region  when  packed  four  per  word. 
The  value  of  IK  is  dependent  on  the  variable 
region  size  as  defined  by  RSS2.  The  order- 
ing of  this  elevation  data  and  its  relation 
to  the  corresponding  grid  locations  within 
the  region  is  the  same  as  that  output  from 
RSS6.  Figure  II. 7. 3 shows  a sample  print- 
out from  RSS7. 

2. 2. 7. 9 Externals . System  routines  called  by  RSS7  are  OPENMS,  READMS, 

WRITMS , CLOSMS,  and  UNIT  which  are  mass  storage  I/O  routines.  RSS7  calls 
no  other  subroutines. 

2.2.7.10  Error  Conditions.  RSS7  has  three  error  messages.  The  first  error 
message  is  output  if  the  resolution  values  of  the  terrain  data  and  planimetry 
data  are  different.  Its  format  is: 

ERROR-CULTURE  RESOLUTION  NOT  EQUAL  TO  TERRAIN  RESOLUTION 
The  second  error  message  is  output  if  an  error  is  found  during  input/output. 

Its  format  is: 


Word  N+l- 
Word  N+UT 


HIT  EOF  OR  PARITY  ON  I/P  TAPE 

The  third  error  message  is  output  if  the  terrain  file  is  wholly  or  partially 
empty.  Its  format  is: 

END  OF  TERRAIN  INPUT  FILE  REACHED  ON  REGION  . 
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2.2.7.11  Program  Flowchart.  Figure  II.  11  shows  the  flow  diagram  of  RSS7. 

2.2.7.12  Program  Listing.  The  complete  listing  of  RSS7  is  attached. 

2.2.8 

2. 2. 8.1  Program  Name . RSS8 

2. 2.8.2  Storage.  Approximately  43K  words  of  core 

2.2.8. 3 Run  Time.  Approximately  1800  cp  seconds  and  45  minutes  wall 
clock  time. 

2. 2. 8. 4 Cost.  Approximately  375  dollars. 

2. 2. 8. 5 Program  Function.  RSS8  creates  the  radar  scan  lines  from  which 
the  final  radar  scene  is  constructed. 

2. 2.8.6  Program  Description.  In  polar  coordinates,  any  point  on  the  map 
may  be  described  by  the  following  expression: 

x = x + r cose 

P t p 

y = + r SIN© 

p t p 

where  (X  . Y ) is  the  location  of  the  target,  r the  radial  distance  from 
t t P 

the  target,  and  0 the  polar  angle.  The  quantity  IFREQ  defines  the  values 
of  0 to  be  used  in  the  generation  of  the  scan  lines.  These  values  are: 

0n  = n/IFREQ,  0^  n £ 360*IFREQ 

and  result  in  the  generation  of  a circle  consisting  of  360*IFREQ  seen  lines 
separated  by  an  angular  increment  of  1/ IFREQ  degrees.  The  radar  range  R 
is  expressed  in  units  represented  by  variable  RESK,  the  grid  spacing  for 
the  map  data.  The  variable  r therefore  has  the  range  O^r^-R.  At  present 
the  spacing  for  terrain  and  planimetry  data  must  be  equal.  However  an 
enhancement  which  allows  for  different  resolution  values  has  been  coded 
into  the  program. 
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I 

With  the  parameters  thusly  defined,  the  generation  of  the  sweep  lines 

proceeds  as  follows.  We  assume  that  tie  initial  value  of  0 =0,  i.e.,  n=0 

n 

and  that  we  begin  each  scan  line  r =1. 

P 

A.  Calculate  COS0  and  SIN0 

n n 

B.  Calculate  X and  Y 

P P 

C.  Calculate  the  region  in  which  (X  , Y ) lies 

P P 

D.  If  the  record  for  this  region  is  already  in  core,  go  to  (G) 

E.  Read  the  record  for  this  region  into  core 

F.  Copy  the  planimetry  strip  information  for  this  region  into 
IK*IK  core  array  which  serves  as  an  image  of  the  region. 

Radar  return  intensity  levels  are  assigned  in  correspondence 
with  the  feature  codes  by  using  an  internal  translation  table. 
The  radar  return  intensity  for  resolution  elements  not  covered 
by  planimetry  data  is  set  equal  to  a predetermined  background 
intensity.  Tbe  result  is  that  array  (I,J)  contains  the  radar 
return  intensity  (from  planimetry  only)  for  the  resolution 
element  at  X=I,  Y=J. 


G.  Calculate  which  resolution  element  contains  (X  , Y ) . Denote 

this  by  ARRAY  (I  , J ) . P P 

P P 

H.  Copy  ARRAY  (I  , J ) and  the  corresponding  terrain  value  from 

P P 

the  input  region  record  into  a linear  array  which  stores  the 
information  for  the  scan  line. 

I.  Increment  r by  1. 

J.  If  riR  go  to  (B) 

K.  Increment  0 by  1/IFREQ 

L.  If  0>36O°  we  are  finished.  If  05360°  go  to  (A) 

The  following  is  a description  of  the  important  variables  in  RLMS8. 

Description 

Starting  angle  of  radar  swap-scan 

Elevation  values  of  each  point  in  the 
terrain  region 


Label 

Type 

ALPHA 

REAL 

ARRAY 

INTEGER 

ARRAY 
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I 

I 

I 

I 

I 

I 


jL 
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Label 

Type 

Description 

ALTTDE 

REAL 

Altitude  of  radar  in  feet  above  ground 
level . 

CONV 

REAL 

Conversion  factor.  Computes  radians 
per  degree 

DCONV, 

ECONV 

REAL 

Conversion  factor.  Computes  grid  points 
per  nautical  mile 

CTO 

REAL 

COS  of  sweep  angle 

DATA 

DLINK 

INTEGER 

ARRAY 

REAL 

ARRAY 

Feature  intensity  level  conversion  table 

Array  for  parameter  input  record  of  REAL 
data  values 

HALFX, 

HALFY 

REAL 

Rounding  factor.  Value  is  .5 

IA,  IAA, 

IB,  IBB 

INTEGER 

Row  and  column  of  current  gross  element 
number  under  beam  scan 

I ANGLE 

INTEGER 

Scan  rotation  angle  counter 

IANX, 
I ANY 

INTEGER 

Polar  position  of  point  within  region 

I ARRAY 

INTEGER 

ARRAY 

Intensity  codes  for  each  point  in  the 
culture  region 

IBACK 

INTEGER 

Intensity  code  for  background  of  picture 

ID 

INTEGER 

Number  of  grid  elements  in  a line  sweep. 
Maximum  value  is  2000. 

I FREQ 

INTEGER 

The  number  of  radial  scan  lines  to  be 
constructed  per  degree.  This  variable  is 
part  of  the  parameter  list. 

IHT 

INTEGER 

Elevation  value  for  point  along  radial 
scan  line. 

II 

INTEGER 

ARRAY 

Array  of  assigned  intensity  codes  for 
each  feature 

IFL 

INTEGER 

Bit  pointer  into  II  array 
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Label 

Description 

IIP 

INTEGER 

Bit  pointer  to  elevation  value 

IIX,  IIY 

INTEGER 

Word  pointers  to  elevation  value 

IK 

INTEGER 

Region  size  for  terrain  data 

IKC 

INTEGER 

Region  size  for  culture  data 

IN 

INTEGER 

Number  of  words  for  region  row  data 

ILIM 

INTEGER 

Number  of  words  in  current  region  record 

IMIN , 
IMAX 

INTEGER 

Beginning  and  ending  scan  angle  (+1) 

INR 

INTEGER 

Equivalent  to  REG  used  as  an  index  for 
READMS. 

IN2 

INTEGER 

Intensity  code  for  point  along  radial 
scan  line 

IN2W 

INTEGER 

Input  buffer  for  two  word  file  output 
by  RSS2. 

IXL,  IYL 

INTEGER 

Row  and  column  element  values 

I TERR 

INTEGER 

Switch  for  terrain  data.  1=  terrain  data. 
0=  no  terrain  data.  This  variable  is  part 
of  the  parameter  list. 

JLIM 

INTEGER 

Number  of  data  values  region  row. 

KX,  KY 

INTEGER 

Row  and  column  of  current  element  under 
beam  scan  within  a given  culture  region 

KKX,  KKY 

INTEGER 

Row  and  column  of  current  element  under 
beam  scan  within  a given  terrain  region 

MSIX 

INTEGER 

ARRAY 

Index  directories  used  by  mass  storage  I/O. 

NLINK 

INTEGER 

ARRAY 

Array  for  parameter  input  record  of  INTEGER 
data  values 

NRX 

INTEGER 

Maximum  number  of  allowable  regions  along 
X-AXIS . 
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Label 

Description 

NREG 

INTEGER 

Region  number  for  which  data  is  currently 
in  core. 

P 

REAL 

Current  angle  of  sweep  line  in  degrees. 

RANGE 

REAL 

Radar  radius  of  coverage  in  nautical  miles. 
This  variable  is  part  of  parameter  list. 

REG 

INTEGER 

Region  number  of  current  sweep  line  element. 

RESK 

REAL 

Resolution  unit  size  in  feet. 

SCNLNE 

INTEGER 

ARRAY 

Scan  line  storage  for  one  row  of  culture 
intensities  and  corresponding  elevations 

ST« 

REAL 

SIN  of  sweep  angle 

THETA 

REAL 

Current  angle  of  sweep  line  in  radians 

WORD 

INTEGER 

ARRAY 

Array  for  two  word  file  output  from  RLMS2 

XI,  YI 

INTEGER 

Radar  location  in  terms  of  gross  data  grid 

X,  Y 

INTEGER 

Target  location 

2. 2. 8. 7 Input.  rss8  requires  three  input  disk  files.  The  files  required 
are  the  two  word  output  file  (Tape  20)  from  RSS2 , the  output  parameter  file 
(Tape  6)  from  RSS7  and  the  random-access  map  data  file  (Tape  12)  from  PS S7, 

2. 2.8.8  Output . RSS8  outputs  one  disk  file  (Tape  4).  This  file  consists 

of  360*IFREQ  records  of  4000  words  each.  The  file  is  sequential  in  nature 

with  the  records  being  ordered  by  their  corresponding  angle  variable  0,  i.e. 

the  record  for  scan  line  at  0=0°  is  the  first  record  on  the  file  while  that 
o 

for  0=360  is  the  last  record. 

The  first  2000  words  contain  the  intensity-of-return  values  for 
each  r,  O^r^R.  Clearly,  if  R<2000  the  remaining  words  are  blank.  Words 
numbered  2000+1,  ill £2000,  contain  the  corresponding  terrain  elevation 


I 


'I 


I 


values.  For  example,  if  IFREQ=2,  then  the  fourth  record  has  the  following 
significance. 

A.  It  corresponds  to  9=4/IFREQ=2° 

B.  The  Ith  word  and  the  I +2000 ^ word  respectively  contain 

the  strength  of  return  and  the  terrain  elevation  for  the  point: 

X = X + I *COS  2° 

P t 

Y = Y + I*SIN2° 

P t 

Ttoe  printout  from  RSS8  is  shown  in  Figure  II. 8. 2 and  is  self-explanatory. 

2. 2. 8. 9 Externals.  System  routines  called  by  RSS8  are  three  mass  storage 
I/O  function  routines  OPENMS,  READMS,  and  UNIT;  and  three  trigonometric 
function  routines  ATAN2,  COS,  and  SIN.  RSS8  calls  no  other  subroutines. 

2.2.8.10  Error  Conditions.  RSS8  has  three  error  conditions.  Error  condi- 
tion number  one  occurs  when  there  are  more  than  500  culture  strips  for  any 
given  region.  The  format  of  the  error  message  for  this  condition  is  as 
follows : 

VALUE  COUNT= , IN  ERROR  FOR  REGION 

....OCCURRED  DURING  SWEEP  ANGLE  OF  DEGREES 

Error  condition  number  two  occurs  when  the  index  for  the  region  data  extends 
beyond  the  region  bounds.  The  error  messages  for  this  condition  are  as 
follows: 

COL  ELEMENT  VALUE= , IN  ERROR  FOR  REGION 

OCCURRED  DURING  SWEEP  ANGI£  OF  DEGREES 

OR 

ROW  ELEMENT  VALUE* , IN  ERROR  FOR  REGION 

....OCCURRED  DURING  SWEEP  ANGLE  OF  DEGREES 


I 
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FIGURE  II  • 1^  - PROGRAM  RSS8  FLOWCHART 
(Page  2 of  2) 


The  third  error  condition  occurs  if  there  is  an  error  found  during  input/ 
output.  The  format  of  its  error  messages  is  as  follows: 

EOF  OR  PARITY  ERROR  ON  UNIT  6,  RSS8 

....OCCURRED  DURING  SWEEP  ANGLE  OF  DEGREES 

OR 

EOF  OR  PARITY  ERROR  ON  O/P  FILE  4 

....OCCURRED  DURING  SWEEP  ANGLE  OF  DEGREES 

OR 

INPUT  ERROR  ON  UNIT  20 

2.2.8.11  Program  Flowchart.  Figure  II.  12  shows  the  flow  diagram  of 
RSS8. 

2.2.8.12  Program  Listing.  The  complete  listing  of  RSS8  is  attached. 

2.2.9 

2. 2.9.1  Program  Name.  RSS9 

2 • 2 . 9 . 2 Storage.  Approximately  10K  words  of  core 

2. 2. 9. 3 Run  Time.  Approximately  1000  cp  seconds  and  4U  minutes  wall 
clock  time 

2. 2. 9. 4 Cost.  Approximately  160  dollars 

2- 2.9. 5 Program  Function.  rss9  applies  the  radar  effects  to  the  radial 
scan  line  data  from  RSS8. 

2. 2. 9. 6 Program  Description.  The  three  radar  effects  considered  by  RSS9  are: 

(1)  Lambert's  Law  effect  which  determines  the  percentage  of  the 
incident  radar  signal  reflected  back  to  the  source. 

(2)  Shadow  effect  which  takes  into  account  the  fact  that  certain 
areas  on  the  ground  may  be  blocked  from  view  by  tall  objects 
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(i.e.  mountain  peaks)  lying  in  the  line-of-sight  from  the 
radar  location  to  the  area  in  question. 

(3)  Altitude  layover  effect  which  takes  into  account  that  the 
radar  perceives  the  location  of  high  altitude  terrain  to  be 
closer  to  the  radar  than  it  actually  is,  and  below  sea  level 
terrain  to  be  farther  from  the  radar  than  it  actually  is. 

RLMS9  contains  a subroutine  used  to  scale  the  final  image  so  that  scenes 
generated  from  different  altitudes  will  all  be  the  same  size  when  displayed 
on  the  DICOMED  plotter.  RSS9  also  contains  a subroutine  which  converts  the 
radial  formatted  points  back  to  cartesian  coordinates  for  eventual  display. 
The  following  is  a description  of  the  important  variables  in  RSS9. 


Label 

Type 

Description 

ALT 

INTEGER 

Altitude  of  the  radar  above  sea  level 
(in  feet) 

ATARG 

REAL 

Altitude  of  the  radar  above  the  target 
(in  feet) 

AVINS 

REAL 

ARRAY 

Array  which  contains  the  average  of 
intensity  for  each  point  along  the 
scan  line,  i.e.,  each  value  in  the  array 
equals where  n is  the  number  of 

points  along  the  scan  line,  and  I is  the 
intensity  value  for  the  point  r under  con 
sideration;  and  360  *IFREQ*720. 

CONV 

REAL 

Conversion  factor.  Computes  radians  per 
degree . 

DATA 

INTEGER 

ARRAY 

Array  for  input  data  record  of  combined 
culture  and  terrain  values  by  sweep  line 

ELEV 

INTEGER 

ARRAY 

The  2001^-4000 th  word  of  the  input  data 
which  contains  the  terrain  values. 

ENTS 

REAL 

The  intrinsic  intensity  value  for  a point 
with  Lambert's  Law  included  as  a factor. 
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Label 

Type 

Description 

I CODE 

INTEGER 

ARRAY 

Array  of  counters.  Counts  the  number  of 
each  color  code  generated  before  Lambert's 
Law  is  considered. 

I COLOR 

INTEGER 

ARRAY 

Array  of  counters.  Counts  the  number  of 
each  color  code  generated  before  altitude 
layover  is  considered. 

IC0L02 

INTEGER 

ARRAY 

Array  of  counters.  Counts  the  number  of 
each  color  code  generated  after  altitude 
layover  is  considered. 

IDIST 

INTEGER 

Number  of  points  along  the  radial  scan  line. 

I FREQ 

INTEGER 

The  number  of  radial  scan  lines  to  be  con- 
structed per  degree. 

IMIN, 
I MAX 

INTEGER 

Beginning  and  ending  scan  angle  (+1) 

INS 

INTEGER 

ARRAY 

The  lSt-2000t^1  word  of  the  input  data  which 
contains  the  culture  values. 

IN2W 

INTEGER 

Input  buffer  for  two  word  file  output  by 
RSS2. 

ISHADE 

INTEGER 

Counter  for  the  number  of  points  computed 
to  be  a shadow  point  by  subroutine  LAMBERT 

I TEMP 

INTEGER 

Save  variable  for  IDIST 

KOUNT 

INTEGER 

Counter  for  the  number  of  radial  sweep 
points  processed. 

MM 

INTEGER 

Sweep  angle  (+1) 

RESK 

REAL 

Resolution  unit  size  in  feet  for  the 
culture  data 

REST 

REAL 

Resolution  unit  size  in  feet  for  the 
terrain  data 

RET 

REAL 

ARRAY 

Array  containing  the  predetermined  inten- 
sity strengths  for  intensity  codes  0-63 

SCALE 

REAL 

Scaling  factor.  Scales  the  radar  picture 
up  or  down.  SCALE=IDIST/1000 
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Label 

Type 

Description 

SR 

REAL 

ARRAY 

Array  for  the  strengths  of  return 
computed  by  subroutine  LAMBERT 

as 

THETA 

REAL 

Current  angle  of  sweep  line  in  radians 

WORD 

INTEGER 

ARRAY 

Array  for  two  word  file  output  by 

RSS2. 

2. 2.9. 7 I nput . RSS9  requires  three  input  disk  files  and  one  input  card 
file.  The  files  required  are  the  two  word  output  file  (Tape  20)  from  RSS2, 
the  output  parameter  file  (Tape  6)  from  RSS7,  and  the  radial  scan  line  data 
file  (Tape  4)  from  RSS8.  The  card  file  is  a one  word  parameter  file  which 
contains  the  altitude  of  the  radar  above  the  target. 

2. 2.9.8  Output . RSS9  outputs  four  disk  files  (Tape  3,  Tape  10,  Tape  12, 
and  Tape  14) . These  files  contain  the  data  for  unsorted  raster  converted 
points  output  by  subroutine  REFRMT.  Each  output  file  contains  180000  two 
word  records.  These  180000  records  represent  data  for  one  quarter  of  the 
picture,  i.e.,  90*IFREQ*1000=180000 . Each  record  represents  data  for  a 
single  point  on  the  image.  Its  format  is  as  follows: 

Word  1 MOD (X,  Y)  coordinate  in  units  represented 

by  RESK 

Word  2 Intensity  values  from  0-63.  Here  intensity 

zero  corresponds  to  a maximum  radar  return 
(white)  while  a 63  corresponds  to  zero  return. 
This  system  is  the  same  as  used  on  the  DICOMED 
plotter. 

The  printout  generated  by  RSS9  is  shown  in  Figure  II.  14  The  input  statis- 
tics indicate  the  distribution  of  color  codes  as  output  by  RSS8.  in  this 
case  the  color  code  2 was  used  for  the  background.  The  output  statistics 
indicated  the  distributions  of  color  codes  output  after  execution  of  subroutine 
LAMBERT  and  then  subroutine  ALTIAV. 
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«LHS9  SUCCESSFUL  lNO. 
INPUT  STATISTICS 

13JU75  PUT  LS  WITH 
b49  PIXtLS  WITH 
li>2  PIXtLS  WITH 
394  P l XI  L S WITH 
3UT“UT  STATISTICS 

■J  PIXLLS  WITH 
lb  PIXLLS  WITH 
29  PIXELS  WITH 
28  PIXELS  WITH 
60  PIXtLS  WITH 
32  PIXELS  WITH 
AS  PIXELS  WITH 
62  PIXCLS  WITH 
73  PIXELS  WITH 
80  PIXELS  WITH 

as  pixels  with 
9b  PIXELS  with 
94  PIXCLS  hITH 
149  PIXELS  WITH 
165  PIXeLS  WITH 
186  PIXtLS  WITH 
239  PIXELS  WITH 
279  PIXCLS  WITH 
31b  PIXELS  WITH 
389  PIXELS  WITH 
S 02  »IXELS  WITH 
558  PIXELS  HITH 
687  PIXELS  HITH 
973  PIXCLS  WITH 
1407  PIXCLS  WITH 
2326  PIXELS  WITH 
3693  PIXELS  WITH 
*♦649  PIXtLS  WITH 
5932  PIXELS  WITH 
3322  PIXELS  WITH 
13391  PIXELS  WITH 
21559  PIXtLS  WITH 
30186  PIXELS  WITH 
15275  »IXELS  WITH 
7U29  PIXELS  WITH 
3836  3 IXFL S WITH 
E418  PIXELS  WITH 
lc32  PIXELS  WITH 
1029  PIXtLS  HITH 
79U  PIXELS  WITH 
594  PIXELS  WITH 
437  PIXELS  WITH 
359  PIXELS  WITH 
331  PIXELS  WITH 
308  PIXELS  WITH 
569  PIXELS  WITH 
383  PIXCLS  WITH 
207  PIXELS  WITH 
235  PIXtLS  WITH 
161  PIXELS  WITH 
14>  PIXELS  WITH 
102  PIXtLS  WITH 

105  » IXFLS  WITH 

106  PIXELS  WITH 
441  PIXELS  WITH 
970  PIXELS  WITH 
924  PIXELS  WERE 


DISPLAY  FILE  nOPPll.lL 
COLOR  COOt.  ? 

color  cnor  il 

COLOR  Cujf  \i 
COLUR  C 1UE  1-, 

COLOR  CTJC  ft 
COLOR  COOt  9 
LOL  Op  CODE  10 
COLOR  COJt  11 
COLOR  COut  12 
COLOP  COlir  13 
COLOR  COUt  14 
COLOR  COJt  15 
COLOR  COOF  lb 
COLOR  C03E  17 
COLOR  CODE  18 
COLOR  CODE  19 
COLOR  CODE  20 
COLOR  COUE  21 
COLOR  COOt  22 
COLOR  CODE  23 
COLOR  COUE  24 
COLOR  CODE  25 
COLOR  CODE  26 
COLOR  CODE  27 
COLOR  COOF  28 
COLOR  COOE  29 
COLOR  COOF  30 
COLOR  COOF  31 
COLOR  COOE  32 
COLOR  COOF  33 
COLUP  COOF  34 
COLOR  COOE  35 
COLOR  COOE  36 
COLOR  COOt  37 
COLOR  COOE  38 
COLOR  CODE  39 
COLOR  COOE  40 
COLOR  COOE  41 
COLOR  COOE  42 
uOLOP  CD Dc  43 
COLOR  COJt  44 
COLOR  COOE  45 
COLOP  COOt  46 
COLOR  CODr  47 
COLOR  COOE  48 
COLOR  COOE  49 
COLOR  COUE  50 
COLOR  COOE  SI 
COLOR  CDOE  62 
COLOR  COOE  53 
COLOR  COOE  54 
COLO®  COJE  55 
COLOR  COOE  56 
COLOR  COOE  57 
COLOR  COOF  56 
COLOR  COOE  59 
COLOR  COOF  69 
COLOR  COOF  61 
COLOR  COOF  o2 
COLOR  COOE  63 
8SSTGNLJ  TO  TRE  SHAOOW 
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FIGURE  11.14  - SAMPLE  PRINTOUT  FROM  RSS9 


FIGURE  II . 15-  PROGRAM  RSS9  FLOWCHART 
(Page  2 of  2) 
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2. 2. 9. 9  Externals . System  routines  called  by  RSS9  are  ALOGIO,  a logarithmic 


function  routine;  ATAN2,  a trigonometric  function  routine;  and  UNIT,  a mass 
storage  I/O  function  routine.  Subroutines  called  by  RSS9  are  LAMBERT,  SIZE, 
REFRMT,  and  ALTLAV.  Each  of  these  subroutines  are  discussed  under  sections 
2.2.9.13  through  2.2.9.16. 

2.2.9.10  Error  Conditions . RSS9  has  two  error  conditions.  The  first 
error  condition  occurs  if  there  is  an  error  during  input  or  output  to  the 
disk  files.  One  of  the  following  messages  is  output  for  this  error  condition: 

RSS9  - ERROR  OR  EOF  ON  UNIT  6 
OR 

RSS9  - INPUT  ERROR  ON  UNIT  20 
or 

RSS9  - EOF  OR  PARITY  ERROR  IN  I/P  FILE  4 
The  second  error  condition  occurs  if  the  resolution  value  for  the  terrain 
data  is  different  from  that  for  the  planimetry  data.  The  message  output 
for  this  condition  is  as  follows: 

ERROR-CULTURE  RESOLUTION  NOT  EQUAL  TO  TERRAIN  RESOLUTION 

2.2.9.11  Program  Flowchart.  The  flow  diagram  for  RSS9  is  illustrated 
in  Figure  11.15. 

2.2.9.12  Program  Listing.  The  program  listing  for  RSS9  is  attached. 

2.2.9.13  Subroutine  Description 

2.2.9.13.1  Subroutine  Name  LAMBERT 

2.2.9.13.2  Sunmary.  LAMBERT  computes  the  Lambert's  Law  and  shadow  radar 


effects . 


2.2.9.13.3  Description  of  Processing.  The  calling  sequence  for  LAMBERT  is 

LAMBERT  (ELEV,  ALT,  SR,  REST) . For  a description  of  the  calling  arguments 

see  the  list  of  important  variables  for  this  subroutine.  The  Lambert's 

Law  effect  is  illustrated  in  Figure  11.16. 

The  mathematical  statement  of  this  law  is  that  the  return  from  any 

given  point  on  the  ground  is  reduced  by  a factor  COS  8 where  8 is  the 

L L 

angle  between  the  incident  radar  beam  and  the  normal  to  the  terrain  at  the 

point  of  interest.  Ihis  is  to  say,  if  I is  the  intrinsic  intensity  of  the 

background  or  of  the  planimetry  feature  located  at  a given  point,  then  the 

actual  radar  return  intensity  is  given  by: 

I = I *COS8 
r L 

The  implementation  of  this  effect  is  quite  straightforward.  For 

each  value  of  r,  the  slope  of  the  line  through  the  points  r-1  and  r+1  is 

calculated  from  the  corresponding  terrain  elevation  values  thus  yielding 

the  angle  8t>  The  vector  perpendicular  to  this  line  defines  the  normal  to 

the  surface  at  r.  The  declination  angle  of  the  radar  beam  8„  can  be  cal- 

D 

culated  from  r,  the  terrain  elevation  at  r,  and  the  radar  altitude.  Using 

these  angles,  the  angle  6 can  be  determined  from  elementary  trigonometry. 

Li 

The  shadowing  effect  is  illustrated  in  Figure  II.  17.  The  nature 
of  the  effect  is  simply  that  points  on  the  ground  may  be  in  the  shadow 
terrain  peaks  lying  between  the  radar  and  the  ground  point  of  interest. 

Such  points  will  be  invisible  to  the  radar  and  must  therefore  be  assigned 
a zero  intensity-of-return. 
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FIGURE  11.1' i-  ILLUSTRATION  OF  SHADOW  EFFECT 


The  following  is  a description  of  the  important  variables  in  LAMBERT. 


Label 

Tjree 

Description 

A 

REAL 

The  square  of  the  declination  angle  of 
the  radar  beam. 

ALT 

INTEGER 

Subroutine  argument.  The  altitude  of 
the  radar  above  sea  level. 

AR 

REAL 

The  minimum  elevation  required  for  a 
ground  point  to  be  no  longer  in  shadow. 

B 

REAL 

The  square  of  the  slope  of  the  line  through 
r-1  and  r+1  of  any  given  point  r. 

ELEV 

INTEGER 

ARRAY 

Subroutine  argument.  Array  of  elevation 
values  for  each  point  along  a radial  line. 

IDIST 

INTEGER 

Number  of  points  per  radial  scan  line 

IR 

INTEGER 

Equals  variable  'R'  with  a cutoff  value 
defined  by  variable  'RD'. 

IS  HADE 

INTEGER 

Number  of  pixels  assigned  to  shadow 

MT 

REAL 

The  slope  of  the  line  through  r-1  and 
r+1  of  any  given  point  r. 

MV 

REAL 

The  declination  angle  of  the  radar  beam 

R 

INTEGER 

Current  point  along  a radial  scan  line 

RD 

INTEGER 

Minimum  value  of  R in  order  to  avoid 
ending  up  with  a bright  area  in  the 
center  of  the  final  image 

RES 

REAL 

Resolution  element  size  in  feet  scaled 
to  fit  new  image  size 

REST 

REAL 

Subroutine  argument.  Resolution  element 
size  in  feet. 

RP 

REAL 

The  calculated  value  of  'R'  if  point  R 
is  outside  the  view  range  of  the  radar 
at  location  r*0. 

SCALE 

REAL 

Scale  factor  for  adjusting  the  image  size 

SHADOW 

REAL 

The  strength  of  return  for  a shadow  point 
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Label 

Type 

Description 

SLOPE 

REAL 

Same  as  'MV' 

SR 

REAL 

Subroutine  argument. 

Strengths 

of 

ARRAY 

return  for  each  radial  point. 

SRT 

REAL 

Computed  strength  of 

return  for 

a point 

before  being  multiplied  by  the  factor 
RD/5C0  to  avoid  getting  a bright  area 
in  the  center  of  the  final  image. 

X REAL  The  factor  RD/500  used  to  avoid  getting 

a bright  area  in  the  center  of  the  final  image 

2.2.9.13.4  Error  Conditions.  None 

2.2.9.13.5  Subroutine  Flowchart.  The  flow  diagram  for  LAMBERT  is  illustrated 
in  Figure  11.18. 

2.2.9.13.6  Subroutine  Listing.  The  listing  for  LAMBERT  is  included  with 
that  of  RLMS9. 

2.2.9.14 

2.2.9.14.1  Subroutine  Name  ALTLAV 

2.2.9.14.2  Summary.  ALTLAV  computes  the  altitude  layover  radar  effect. 

2.2.9.14.3  Description  of  Processing.  The  calling  sequence  for  ALTLAV  is 
ALTLAV (ELEV,  ALT,  INS,  ATARG,  RET,  SR,  REST) . For  a description  of  the 
calling  arguments  see  the  list  of  inqaortant  variables  for  this  subroutine. 

The  altitude  layover  effect  shifts  the  location  of  the  terrain. 

Ihe  length  and  direction  of  the  shifting  of  each  processed  point  of  the 
terrain  depends  on  the  height  of  that  point.  This  effect  is  implemented 
by  considering  each  point  R along  a radial  line  and  computing  a new  altitude 
layover  point  Rl.  Given  (1)  the  distance  of  the  point  (R)  away  from  the 
target  point  along  a radial  line,  (2)  the  height  of  the  radar  above  sea 
level  (A),  and  (3)  the  elevation  of  the  point  R above  sea  level  (E) ; the 


2-81 


distance  (D)  from  the  radar  to  the  elevated  point  along  the  radial  can 
be  computed  as  follows: 

i — 2 — r 

D = y (A-E  + R 

Using  a principle  of  elementary  trigonometry,  we  know  that  a tangent  line 
drawn  the  length  of  the  elevation  of  point  R to  some  point  R1  at  target 
level  along  the  radial  line  yields  a distance  D1  equal  in  length  to  distance 
D.  Hie  point  R1  can  be  computed  similarly  to  the  above  equation.  Having 
computed  D=D1,  and  given  the  radar  altitude  ablve  target  level  (ALT),  the 

point  Rl  computed  at  target  level  is  as  follows: 

2 2 2~ 

ALT  + (A-E)  + R 


The  following  is  a description  of  the  important  variables  in  ALTLAV: 


Label 

Tyge 

Description 

ALT 

INTEGER 

Subroutine  argument.  The  altitude  of 
the  radar  above  sea  level. 

ATARG 

REAL 

Subroutine  argument.  The  altitude  of 
the  radar  above  the  target 

OX 

REAL 

The  difference  in  feet  between  the  dis- 
tance of  a reference  point  from  the 
target  location  and  the  calculated 
altitude  layover  distance  from  the 
target. 

EI£V 

INTEGER 

ARRAY 

Subroutine  argument.  Array  of  elevation 
values  for  each  point  along  a radial  line 

ENTS1 

REAL 

The  density  value  for  color  code  0. 

ENTS 

REAL 

ARRAY 

Array  containing  the  density  values  for 
each  point  along  radial  line.  When  alti- 

tude  layover  requires  two  density  values 
to  be  added,  the  sum  is  placed  into  this 
array . 
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Label 

Type 

Description 

IDIST 

INTEGER 

Number  of  points  per  radial  scan  line 

INS 

INTEGER 

ARRAY 

Subroutine  argument.  Intensity  codes  for 
each  point  along  a radial  scan  line 

I SHADE 

INTEGER 

Number  of  pixels  assigned  to  shadow 

RES 

REAL 

Resolution  element  size  in  feet  scaled 
to  fit  new  image  size 

REST 

REAL 

Subroutine  argument.  Resolution  element 
size  in  feet 

RET 

REAL 

ARRAY 

Subroutine  argument.  Density  values  of 
each  of  the  64  color  codes  used 

RP 

REAL 

The  calculated  altitude  layover  distance 
in  feet  from  target 

RPM 

REAL 

Radial  distance  of  reference  point  in 
feet  from  target 

R1 

INTEGER 

Grid  location  along  radial  of  the  calcu- 
lated altitude  layover  point 

R2 

INTEGER 

Maximum  value  of  r for  points  defined 
by  (r,  9)  that  a radar  positioned  at  r=0 
can  view 

R3 

INTEGER 

Minimum  value  of  r for  points  defined  by 
(r,  9)  that  a radar  positioned  at  r=0 
can  view 

SR 

REAL 

ARRAY 

Subroutine  argument.  Strengths  of  return 
for  each  point  calculated  by  subroutine 
LAMBERT 

SCALE 

REAL 

Subroutine  argument.  Scale  factor  for 
adjusting  the  image  size 

XP 

REAL 

The  distance  in  feet  of  a radial  point 
from  the  vertex  minus  the  value  of  DX 

X 

REAL 

The  calculated  altitude  layover  distance 
in  feet  from  the  target  location 
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FIGURE  11.19  - SUBROUTINE  ALTIAV  FLOWCHART 
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2.2.9.14.4  Error  Conditions.  There  is  only  one  error  condition  which  will 
cause  an  error  message  output.  That  condition  is  when  the  color  code  is 
computed  to  be  less  than  zero.  When  this  occurs,  the  erroneous  code  is 
printed  and  processing  continues. 

2.2.9.14.5  Subroutine  Flowchart.  The  flow  diagram  for  ALTLAV  is  illustrated 
in  Figure  11.19. 

2.2.9.14.6  Subroutine  Listing.  The  listing  for  ALTLAV  is  included  with 
that  of  RSS9. 

2.2.9.15 

2.2.9.15.1  Subroutine  Name  SIZE 

2.2.9.15.2  Summary . SIZE  scales  the  final  image  up  or  down  so  that  scenes 
generated  from  different  altitudes  will  all  be  the  same  size  when  displayed 
on  the  DICOMED  plotter. 

2.2.9.15.3  Description  of  Processing.  The  calling  sequence  for  SIZE  is 
SIZE  (ELEV,  INS) . For  a description  of  the  calling  arguments  see  the  list 
of  important  variables  for  this  subroutine.  The  need  for  this  operation  is 
illustrated  in  Figure  II.  20,  and  is^due  to  the  fact  that  with  a fixed  radar 
aperture  the  ground  distance  covered  by  the  radar  will  be  the  function  of 
the  altitude  of  the  radar.  Since  we  will  associate  one  resolution  element 
with  one  pixel  on  the  DICOMED  display,  a smaller  ground  range  will  result 

in  a progressively  smaller  display.  Subroutine  SIZE  therefore  confutes  a 
scale  factor  equal  to: 

F = 1000/NO.  OF  RESOLUTION  ELEMENTS  IN  THE  GROUND  RANGE 
We  have  chosen  1000  resolution  elements  as  the  desired  DICOMED  image  radius. 
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covering  a ground  range  of  19.194  nautical  miles  for  a radar  located  32000 
feet  above  ground  level.  If  F is  computed  to  be  negative  the  image  is 
scaled  up;  if  positive,  the  image  is  scaled  down.  This  is  done  simply  by 
assigning  the  elevation  and  radar  return  for  distance  r to  distance  F*r. 

For  points  between  F*r  and  F*(r+1),  return  intensities  are  set  equal  to  the 
background  value  and  elevations  are  computed  by  linear  interpolation. 

The  following  is  a description  of  the  important  variables  in  SIZE. 


Label 

Type 

Description 

ELESV 

INTEGER 

ARRAY 

Holding  area  for  elevation  values  that  are 
being  repositioned  and  scaled  to  new  image 
size 

ELEV 

INTEGER 

ARRAY 

Subroutine  argument.  Array  of  elevation 
values  for  each  point  along  a radial 
scan  line 

IDIST 

INTEGER 

Number  of  points  per  radial  scan  line 

IE 

INTEGER 

Pointer  to  where  to  place  scaled  image 
data  in  the  holding  areas 

INS 

INTEGER 

ARRAY 

Subroutine  argument.  Array  of  intensity 
codes  for  each  point  along  a radial  scan 
line 

INSV 

INTEGER 

ARRAY 

Holding  area  for  intensity  values  that 
are  repositioned  for  the  new  image  size 

ISHADE 

INTEGER 

Number  of  pixels  assigned  to  the  shadow 

J 

INTEGER 

Pointer  to  elevation  value  to  be  scaled 
for  new  image 

KK 

INTEGER 

Pointer  to  intensity  code  to  be  retrieved 
for  new  image 

SCALE 

REAL 

Scale  factor  for  adjusting  the  image  size 

FIGURE  11.21-  SUBROUTINE  SIZE  FLOWCHART 
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2.2.9.15.4  Error  Conditions.  None 

2.2.9.15.5  Subroutine  Flowchart.  The  flow  diagram  for  SIZE  is  illustrated 
in  Figure  11.21. 

2.2.9.15.6  Subroutine  Listing.  Hie  listing  for  SIZE  is  included  with 
that  of  RSS9. 

2.2.9.16 

2.2.9.16.1  Subroutine  Name.  REFRMT 

2.2.9.16.2  Summary . REFRMT  converts  the  radial  data  back  to  raster 
cartesian  coordinates. 

2.2.9.16.3  Description  of  Processing.  The  calling  sequence  for  REFRMT 
is  REFRMT ( DATA , THETA).  For  a description  of  the  calling  arguments  see 
the  list  of  important  variables  for  this  subroutine.  For  purpose  of  dis- 
play on  the  DICOMED  plotter,  the  image  information  must  be  in  raster  for- 
mat. REFRMT  performs  the  preliminary  conversion  from  radial  to  raster. 
Processing  one  radial  line  at  a time  REFRMT  converts  each  radial  point 

to  cartesian  coordinates  as  follows : 

YX  = 2000* ( (Y-l) +X) 

This  YX  value  with  its  corresponding  color  code  is  output  as  a record  for 
a single  point  into  a disc  output  file.  When  all  the  points  have  been 
processed,  the  result  will  be  four  disk  files,  each  of  which  contains  1/4 
of  the  image.  The  following  is  a list  of  the  important  variables  of 
this  subroutine. 


Label 

1 

Description 

CTH 

REAL 

Cosine  of  the  radial 

angle. 

DATA 

INTEGER 

Subroutine  argument. 

Array  containing 

ARRAY 

elevations  and  intensity  codes. 
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Label 

Type 

Description 

9 

ID 

INTEGER 

Current  point  along  radial  line. 

D 

REAL 

ff  H f*  It  II 

IDIST 

INTEGER 

Number  of  points  per  radial  scan  line. 

I LINE 

INTEGER 

Number  of  points  per  raster  line. 

ISHADE 

INTEGER 

Number  of  pixels  assigned  to  shadow. 

IUNIT 

INTEGER 

Output  device  unit  number. 

KOUNT 

INTEGER 

Counter  of  number  of  points  processed 
by  this  subroutine. 

NMAX 

INTEGER 

Maximum  allowable  raster  points  per  line. 

NPIS 

INTER 

Same  as  IDIST. 

PI 

REAL 

Mathematical  value. 
= 3.1515928 

! 

OUT 

INTEGER 

ARRAY 

Output  buffer  for  raster  converted 
records . 

| 

STH 

REAL 

Sine  of  the  radial  angle 

| 

SCALE 

REAL 

Scale  factor  for  adjusting  the  image 
size. 

I 

THETA 

REAL 

Subroutine  argument.  Angle  of  current 
radial  line. 

1 

Y 

REAL 

Random  number  between  0-1  generated  by 
the  random  function  routine  RANF. 

1 

2.2.9.16.4  Error 

Conditions.  There  is  only  one  error  condition.  This 

| 

error  condition  occurs  when  there 

is  an  error  found  while  outputting  a 

record  to  the 

disk  file. 

j 

2.2.9.16.5  Subroutine  Flowchart. 

The  flow  diagram  for  REFRMT  is  illustra- 

) 

ted  in  Figure 

II. 

22. 

I 

2.2.9.16.6  Subroutine  Listing. 

The  listing  for  REFRMT  is  included  with 

1 

! 

that  of  RSS9. 
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2.2.10.1  Program  Name.  SORT 

2.2.10.2  Storage.  Approximately  130K  words  of  core. 

2.2.10.3  Run  Time.  Approximately  1750  cp  seconds  and  2 hours 
wall  clock  time. 

2.2.10.4  Cost.  Approximately  350  dollars. 

2.2.10.5  Program  Function.  Program  SORT  sorts  the  point  records  on  each 
of  the  four  files  output  by  RSS9  and  merges  them  into  a single  file. 

2.2.10.6  Program  Description.  SORT  uses  the  CDC  SORT/MERGE  package  to 
combine  the  four  files,  output  by  subroutine  REFRMT,  into  a single  ordered 
file.  This  file  is  ordered  by  a mod  function  of  Y and  X (MOD(Y,C) ) . The 
procedure  followed  is  to  individually  sort  each  of  the  four  files  and 
then  merge  them  into  one  large  ordered  file. 

2.2.10.7  Input . SORT  requires  the  four  output  disk  files  (TAPE3,  TAPE10, 
TAPE12,  and  TAPE14)  from  RSS9. 

2.2.10.8  Output . SORT  outputs  a single  disk  file  (TAPE27)  which  contains 
1.44  million  ordered  records  describing  each  point  to  later  be  output  to 
the  DICOMED.  The  record  format  is  the  same  as  that  from  RSS9. 

2.2.10.9  Program  Listing.  The  program  listing  for  SORT  is  attached. 

2.2.11 

2.2.11.1  Program  Name.  RSS10 

2.2.11.2  Storage.  Approximately  3K  words  of  core. 

2.2.11.3  Run  Time.  Approximately  510  cp  seconds  and  zo  minutes  wall 
clock  time. 

2.2.11.4  Cost.  Approximately  65  dollars. 

2.2.11.5  Program  Function.  RSS10  uses  the  sorted  data  from  SORT  to 
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format  a plot  tape  for  the  DICOMED, 

2.2.11.6  Program  Description.  The  DICOMED  Screen  consists  of  a 2048  x 
2048  grid  with  the  intensity  of  each  grid  element  being  specified  by  a 
6-bit  color  code.  In  the  mode  of  operation  presently  being  employed, 
plotting  is  done  sequentially  in  horizontal  rows  beginning  at  the  top  of 
the  screen.  Therefore,  in  order  to  describe  a picture  it  is  necessary 
to  format  2048  records  of  205  words  each  e.g.,  one  record  for  each  row 
consisting  of  2048  6-bit  color  codes  packed  10  per  word. 

The  present  image  format  calls  for  the  radar  scene  to  be  displayed 
as  a circle  of  radius  1000  pixels.  Therefore,  the  first  step  is  to  generate 
11  blank  (all  white)  records  to  describe  the  top  margin  of  the  picture. 

The  actual  image  begins  in  row  12. 

For  each  of  the  2000  records  containing  image  information,  those 
pixels  actually  lying  in  the  radar-scene  circle  are  colored  black.  As 
should  be  evident,  not  all  points  within  the  circle  will  contain  image 
information  - only  the  points  on  the  radial  scan  lines  contain  data  and 
if  an  angular  spacing  of  >5°  is  used,  these  comprise  only  about  1/4  of 
the  total  number  of  possible  image  points.  Therefore,  the  choice  of  the 
"fill"  color  is  important  and  black  is  chosen  since  it  will  not  intro- 
duce any  correlation  error  when  the  generated  image  is  compared  to  the 
more  dense  output  of  a real  radar. 

The  information  from  the  input  file  is  then  written  over  the  fill 
color.  Since  the  information  is  sorted,  taking  the  points  in  the  order 
they  appear  permits  the  records  to  be  filled  sequentially  left  to  right, 
with  a change  in  Y signalling  the  end  of  a given  line.  When  the  radar 
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circle  is  completed,  a bottom  margin  of  12  blank  records  is  written  to 
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the  plot  tape. 

This  program  also  places  four  fiducial  marks  located  at  each 
edge  of  the  image.  A small  crossmark  is  also  placed  at  the  center  of 
the  image  to  mark  the  target  location. 

The  following  is  a description  of  the  important  variables  in 

RIMS10. 


Label 

Type 

Description 

BACK6 

INTEGER 

The  background  fill  code. 

DATA 

INTEGER 

ARRAY 

Input  buffer  record  containing  infor- 
mation for  a single  point  to  be  plotted. 

IDELTA 

INTEGER 

Half  the  number  of  points  covered  by 
the  bounds  of  radar  circle  for  this  line. 

IDIST 

INTEGER 

Number  of  points  per  radial  scan  line. 

IF  LAG 

INTEGER 

Control  flag.  Value  of  this  flag  deter- 
mines location  of  program  jump. 

I LINE 

INTEGER 

Current  line  number. 

IMAX 

INTEGER 

Right  boundary  point  of  circle  for  this 
line. 

IMIN 

INTEGER 

Left  boundary  point  of  circle  for  this 
line. 

INDEX 

INTEGER 

Line  index  for  top  and  bottom  margins. 

IPASS 

INTEGER 

Control  flag.  Counts  the  number  of  data 
lines  processed. 

IX 

INTEGER 

X - value  of  current  point. 

IY 

INTEGER 

Y - value  of  current  point. 

JBSHIFT 

INTEGER 

ARRAY 

Table  containing  number  of  bits  to  shift 
6-bit  color  code  into  a word  containing 
10  codes. 
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Label 

Type 

Description 

JBYTE 

INTEGER 

Pointer  to  a 6-bit  byte  in  a 10  byte 
word  of  output  buffer. 

1 

JBYTE1 

INTEGER 

Pointer  to  the  byte  containing  the 

left  boundary  point  data  for  current  line. 

1 

JBYTE2 

INTEGER 

Pointer  to  the  byte  containing  the  right 
boundary  point  data  for  current  line. 

1 

JMASK 

INTEGER 

ARRAY 

Array  of  6-bit  masks  to  insert  each  of 
the  ten  color  codes  into  a word. 

1 

JPIXEL 

INTEGER 

The  color  code  for  the  current  point. 

| 

JWORD 

INTEGER 

Pointer  to  a word  in  the  205  word  out- 
put buffer. 

1 

JWORD1 

INTEGER 

Pointer  to  the  word  containing  the  left 
boundary  point  data  for  current  line. 

1 

JWORD2 

INTEGER 

Pointer  to  the  word  containing  the  right 
boundary  point  data  for  current  line. 

1 

LBUF 

INTEGER 

ARRAY 

Array  of  counters  for  each  of  the  64 
color  codes. 

1 

LINBUF 

INTEGER 

ARRAY 

Output  line  buffer. 

1 

LINE 

INTEGER 

Line  counter. 

■ 

MP 

INTEGER 

Output  file  parity. 

1 

NBDRLN 

INTEGER 

Number  of  lines  in  the  top  and  bottom 
margins. 

1 

NPTS 

INTEGER 

Number  of  points  per  scan  line. 

1 

PICSIZ 

INTEGER 

Number  of  points  per  scan  line. 

■ 

2.2.11.7  Input.  RSS10  has  one  input  disk  file  (TAPE1)  containing  the 
sorted  data  from  the  preceding  SORT/MERGE  program. 

2.2.11.8  Output . RSS2  outputs  the  DICOMED  plot  file  (TAPE3)  on  magnetic 
tape.  Each  205  word  record  of  the  file  contains  color  code  data  for  one 
raster  line  packed  ten  6-bit  codes  per  word.  The  final  image  generated 
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by  the  tape  is  presented  in  Figure  II.  23. 


2.2.11.9  Externals . System  routines  called  by  RSS10  are  the  square 
root  function  (SQRT) ; and  the  mass  storage  I/o  function  routine  (UNIT) . 

There  are  no  subroutine  calls  from  this  program. 

2.2.11.10  Error  Conditions.  There  is  one  error  condition  for  RSS10. 

This  error  condition  occurs  when  there  is  an  error  found  while  reading 
or  writing  to  or  from  disk  storage.  The  messages  generated  for  this 

error  condition  are  as  follows.  ■ 

PARITY  ERROR  WHILE  WRITING  TOP  MARGIN  | 

or 

PARITY  ERROR  ON  INPUT  FILE  j 

or 

OUTPUT  PARITY  ERROR. 

2.2.11.11  Program  Flowchart.  The  flow  diagram  for  RSS10  is  illustrated 
in  Figure  11.24. 

2.2.11.12  Program  Listing.  The  program  listing  for  RSS10  is  attached. 

I 

I 

I 

l 

1 

J 

2-9e  r 


FIGURE  11.24-  PROGRAM  RSS10  FLOWCHART 
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3.  Program  Listings 

The  enclosed  program  listings  represent  the  latest  version  of  the 
RSS  software  and  differ  from  the  preceding  documentation  as  follows: 

Program  RSS8 , in  present  form,  will  skip  every  other  point 
along  the  radial  lines  that  it  generates,  degrading  the  radial  resolution 
by  a factor  of  2.  This  was  done  when  the  azimuthal  resolution  was  in- 
creased from  *j°  to  k°  to  keep  the  total  quantity  of  data  constant.  The 
program  also  prints  out  the  X and  Y errors  that  result  when  the  target 
location  is  rounded  off  to  the  center  of  the  nearest  resolution  element. 

Program  RSS9  now  accepts  up  to  four  input  card  parameters, 
via  a NAMELIST  read,  and  degrades  the  radial  resolution  by  a factor  of 
2 as  per  the  changes  to  RSS8.  The  first  input  parameter  is  the  altitude 
of  the  reference  scene,  REFALT=32000.  The  altitude  layover  effect  will 
be  applied  if  the  input  parameter  LAYOVER=0 . The  layover  effect  will 
not  be  applied  if  LAYOVER  +0.  The  default  value  in  the  program  is 
LAYOVER=0 . The  size  of  the  final  output  raster  image  is  specified  by 
the  parameter  IRASTER,  with  a default  value  of  301  (301  pixels  by  301 
lines) . And,  finally,  the  parameter  ZANGLE  (default  = 0.0)  will  allow 
rotation  of  the  image  by  ZANGLE  degrees  in  the  counterclockwise  direction 
before  conversion  to  raster  format. 

Program  RSS10  inputs  five  card  parameters  with  a NAMELIST 
read.  IRASTER,  with  a default  value  of  301,  specifies  the  input  raster 
image  size.  IBLOWUP,  with  a default  value  of  5,  enlarges  the  raster  image 
IBL0MUP  number  of  times  to  better  match  the  output  plotter's  resolution. 
IOFX  and  IOFY  are  X and  Y offsets  that  will  reposition  the  enlarged  out- 
put raster  by  a maximum  of  + IBLOWUP/2  times,  and  these  are  calculated 
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from  the  X and  Y errors  printed  out  by  RSS8.  The  use  of  IBLOWUP,  IOFX, 
and  IOFY  allows  a low  resolution  but  high  geometric  centering  accuracy 
of  the  final  scene.  The  last  parameter  is  ISKIP,  normally  set  equal  to 
IRASTER,  which  tells  the  program  the  maximum  number  of  skipped  pixels 
along  a final  raster  output  line  to  be  filled  in.  The  radial  line  to 
raster  conversion  will  sometimes  leave  raster  pixels  unfilled,  and  the 
program  will  place  in  these  locations  the  average  of  the  intensities  on 
either  side. 
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